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ABSTRACT 


A new program for calculating dual reflector antenna 
radiation patterns has been developed adding one more option 
to the original program developed jointly by NCSU and NASA. 
The previous program was capable of computing patterns for 
single reflector antennas with either smooth analytic sur** 
faces or with surfaces composed of a number of panels. 

Techniques based on the geometrical optics (GO) approach 
are used in tracing rays over the following regions: 

1) From a feed antenna to the first reflector 
surface (subreflector) . 

2) From this reflector to a larger reflector surface 
(main reflector) . 

3) From the main reflector to a mathematical plane 
(aperture plane) in front of the main reflector. 

The equations of GO are also used to calculate the 
reflected field components for each ray making use of the 
feed radiation pattern and the parameters defining the sur- 
faces of the two reflectors. These resulting fields form 
an aperture distribution which is integrated numerically to 
compute the radiation pattern for a specified set of angles. 

Spillover, diffraction and other factors [2] that affect 
the accuracy of the calculation of the far-out sidelobes, 
are neglected. 
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Examples and all test cases eure mentioned to support the 
validity of the new algorithm* 
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1. INTRODUCTION 


The objective of the work reported herein was to develop 
an algorithm to calculate the radiation patterns of Casse- 
grain antennas, which belong to the general group of dual 
reflector emtennas. (See Appendix A.) The approach taken 
is to adopt and extend an existing algorithm which was 
developed for single reflector antennas. 

The original algorithm for single reflector antennas 
was published in 1976 [l] . Later on that year this program 
appeared as a NCSU report [2] , but in a modified version. 
Between 1976 and 1978 this algorithm was extended to deal 
with new surfaces such as ellipsoids and spheres Q3] . In 
1980, Botula modified the algorithm giving it the capability 
to analyze antennas with either smooth analytic surfaces or 
with surfaces composed of a number of panels [s] . 

The method of the electric vector potential and the 
geometrical optics approach were used to compute the radi- 
ation field of the antenna in question. 

This thesis includes: 

1) All modifications and additions inserted into the 
program to increase the accuracy of the calculated results 
for multipanel single reflector antennas; 

2) The equations written to describe hyperbolic sur- 
faces ; and 

3) The equations used to describe all reflections of 
rays from both surfaces of a Cassegrain antenna and the 


intersections of these rays with tae two surfaces. 

FORTRAN G level was the language used in writing the 
algor ithni. The computing time was slightly increased due 
to the fact that more ray tracing is involved in a dual 
reflector antenna case. 
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2. ANALYSIS AND FORMULATION 
2.1 Theoretical Development 

The majority of operations in this algorithm are 
essentially the same as those in the single reflector 
algorithm. The GO approach is applied to calculate the 
reflected electric field using the feed radiation pattern 
and all parameters defining the surfaces comprising a 
reflector antenna. The electric field is computed over 
a planar aperture in front of the reflector surface. As 
a result, an integration over the aperture plane yields 
the radiation patterns of the antenna in question. 

To understand the line of thought emd development of 
the new algorithm it is necessary to review some aspects of 
the old program and see where the new additions appear. A 
more refined and detailed expl^ulation of all equations in 
the old algorithm is given in references [l] to [S] . 

Figures 2.1 and 2.2 depict the coordinatv? systems used 
in the single and dual reflector algorithms. 

The first difference is that the new algorithm has the 
capeOsility of analyzing both dual reflector antennas and 
single reflector antennas, i.e., the old algorithm became 
part of the new one. The two reflector surfaces are des- 
cribed in terms of the reference coordinate system (x, y, 
and z) in which most of the mathematical operations are 
performed. The second difference between the old 2 md new 
programs lies in the types of reflector surfaces that can 
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be analyzed. Previously/ five types were available t 

planes, spheres, ellipsoids, paraboloids, and parabolic 
\ 

cylinders, whereas now hyperboloids can also be treated 



It should be stressed here that these six types of surfaces 
are available for each reflector for the case of dual re- 
flector antennas. 


MAIN 
REFLECTOR 


FEED, 



APERTURE 
^ PLANE 




REFERENCE 
SYSTEM " 


\ 


subreflector 

Fig. 2.2. Coordinate system for a dual 
reflector antenna system 


Spherical coordinates are used for the radiation pat- 
term calculations. The convention used concerning the 
angles 6 and ^ is shown in Figure 2.3. 


xf 



Fig. 2.3. Convention used for angles 6 and ^ 

The feed position is expressed in terms of the primed 
coordinates th* , y'*, and s'*. The feed radiation pattern is 
expressed in spherical coordinates, based on the feud carte* 
sian coordinate system using the same convention for the 
angles e'' and as the reference spherical system. Here, 

6" and are referred to the feed coordinate system. The 
phase center of the feed antenna is the origin of its coor- 
dinate system. 

The two coordinate systems are related to each other 
via a three-dimensional rotational matrix [A] , whose deri- 
vation can be found in [2] . The rotational operation of 
this matrix is used to make the feed system parallel to the 
reference system, making use of the three angles ALPHA, BETA, 
and GAMMA as shown in Figure 2.4. All counterclockwise rota- 
tions are defined as positive when looking in the negative 
direction along the axis of rotation. 
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ALPHA is the rotation about the z^-axis, BETA is the 
rotation about the x^-eucis euid GAMMA is the rotation eUaout 
the V -axis. 


Each ray starts from the feed and is traced up to the 
aperture pleuie. Five pieces of information are associated 
with each ray; a set of angles 0" and , the appropriate 
Q' ^urld polarized electric field strengths and the ini- 
tial phase, all teUcen from the feed antenna pattern. Fig- 
ures 2.5 and 2.6 show all vector operations involved in ray 


A POINT ON 
APERTURE 
PLANE 

♦ Z 


X 



Fig. 2.5. Vector operations for a single reflector antenna 
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The symbols in these figures are defined as follows: 

/V 

1) 3 ^ is a unit vector in the direction of an arbi- 
trary ray incident on the reflector (or on the 
subreflector) . 

2) R is the distance from the phase center of the feed 
to the point at which the incident ray strikes the 
reflector (or the subreflector) . 

A. 

3) n^ is the unit normal vector to the reflector sur- 
face (or the subreflector) . 

A 

4) Sj. is a vector in the direction of the reflected 
ray, (or reflected from the subreflector) and inci- 
dent on the main reflector in the case of a dual 
reflector emtenna. 

5) RM is the distance from(xQ,yQ ,Z q) on the subreflec- 
tor to[xQ2»yQ2^^02^ main reflector, i.e., the 

distance from a point on the subreflector to a 
point at which the reflected ray strikes the main 
reflector. 

A 

6) s ^2 is ^ vector in the direction of the ray reflect- 
ed by the main reflector. 

7) D is the distance from the point of reflection 

to the aperture plane for a single re- 
flector or from the point <XQ2»yo2'^02^ main 

reflector to the aperture pleme for the dual 


reflector case. 
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The unit vector which is expressed in spherical 
feed coordinates is written in its cartesian coordinate 
system as: 

yv A y\ A 

Si = s; X + Sy y + s; z 

where 

s' = sin Q' cos 

A 

Sy = sin Q' sin and 

s' ~ cos %' 

z 

Q' and are also expressed in terms of the feed cartesian 
coordinates. The feed system is not only rotated but trans- 
lated with respect to the reference system. That means that 
a rotation as well as a translation should be performed to 

A 

express the vector in the reference system. To achieve 

this task, the origin of the reference system must be known 

in the feed system. 

A 

The intersection of a ray having the unit vector s^, 

with the reflector or subreflector surface is defined by a 
vector ^ Figure 2.7. 
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Thus ^ s r 8^ - 0^0 provided that and 0^0 are expressed 
in the reference coordinate system. To accomplish the trans 
formation a 3x2 matrix is formed. This matrix has the 

A 

ray \init vector (s^) and the translation vector as its 
columns. The rotational operation takes place by premulti- 
plying [bb] by the rotation matrix [a] . 

[A] [BB] = [B] 


Each ray is now described in the reference system by the 
parametric equations 

X = B^j_r - Bj_2 


y = B 23 _r - B 


22 


z 


B3ir 


B 


32 


The point of intersection is found by solving simultaneously 
the equations mentioned above and the equation of the reflec 

A 

tor surface. To find a vector (s^) in the direction of the 
reflected ray, the unit normal to the reflector surface, at 
the incident point is evaluated and Snell's Law is used, i.e 


A A 


A A A 


Similarly, the reflected field except for phase, is given 
by 


-► -* 

- 2(nQ*E^) Uq - E^ where E^ is the incident field, 

attenuated, of course, by a factor since we assume that 

the reflector is in the far field of the feed antenna. 

All vector operations are the same for both the single 

and dual reflector antenna options. 
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The two options are now considered separately. 

A) Dual Reflector System 

A 

The parametric equations for a ray along s^, which is 
treated now as the incident ray on the main reflector, are: 

X ■ Xq + h cosax 

y = y« + h cosay 


z = Zq + h cosaz 

where h is the distance travelled from the point (xQ,yQ,ZQ) 
on the subreflector along the ray, and 

s 


cosax = 


rx 

S'*’ 

r 


cosay 


= !£Z- 


S, 


cosaz - 


rz 

®r 


direction cosines 


emd s , s , s cure the components of the reflected vector 
rx ry rz ^ 

s^. To find the intersections of the ray and the main re- 
flector, simultaneous solution of the above peurametric 
equations with the equations of the surface of the main 
reflector is required. 

The unit normal to the surface is evaluated at this 
point and used to compute a vector in the direction of the 
reflected ray, i.e., 

®r2 * ®i2 ■ ^ ^”02’®i2^ *^02 

where s ^2 = is a unit vector incident on the main reflector. 
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and Hq 2 Is the unit normal on the surface of the main re- 
flector in ceurtesian components. 

®r2 “ * [^ix2 "^“x02^“x02*^ix2 *^y02 *^iy2 *^202 *®iz2!] 

y [®iy2 "^“y02^"x02*^ix2 ”y02 *^iy2 ”z02 *®iz2|] 

® jj®iz2 “^^z02^”x02*^ix2 \02 *^iy2 - *'z02 *®iz2^] 

where 

®ix2 ^ ®rx 

s . ^ * s 
iy2 ry 

®iz2 ' ®rz 

are the components of the ray vector reflected by the sub- 
reflector. Now if 


®rx2 ' 

®ix2 ”^''x02 

^”x02 

•^ix2 

o 

to 

•®iy2 

+ 

*^z02 

*®iz2 

®ry2 “ 

®iv2 ■^*'y02 

^’^x02 

•®ix2 * 

CM 

O 

>1 

•®iy2 

+ 

^z02 

•®iz2 

®rz2 ^ 

®iz2 "^^z02 

^"x02 

•®ix2 

^y02 

•®iy2 

+ 

"z02 

•®iz2 

then 




A 
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Similarly, the reflected field (Figure 2.8), assuming a per- 
fectly conducting reflector, is given by: 


®r2 “ ^ *"02'®i2’ “02 ■ ®i2 


where = jg . 


E ^2 incident electric field on the main reflec- 

-► 

tor and is the electric field reflected by the subreflec- 
tor. It is seen here that E^ is multiplied by a factor 1/RM 
since the main reflector is assumed to be in the far field 
of the subreflector. 

In component form, 


E 



+ y 


- E 

_£Z 

RM 



= X y ^ ®iz2 

and E^2 becomes 

®r2 ^ ^ [^*'x02^”x02*^ix2'^^y02‘®iy2'^*^z02’^iz2^"®ix2] 

y [^^y02^’^x02’^ix2''’\02*®iy2‘''*'z02’®iz2^"®iy2] 

^ [^*'z02^^x02’^ix2'^\02*®iy2'^”z02*^iz2^"®iz2] 

The procedure of finding the intersection of the reflected 
ray (by the main reflector) and the aperture plane is as 
follows: 

Find the parametric equation for a line along s^2 


given by: 
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X = Xq 2 + h'cosa'x 


y “ ^02 **■ 


where 


z + Zq 2 + h''cosa"z 


cosa X = 


* ’ ®r2 ®rx2 


*r2 


’r2 


y . s 




cosa y = — T 


r2 


s 


r2 


cosa z = 


^ • ®r2 ®rz2 


’r2 


’r2 


and h' is the distance travelled along the ray. The aper- 

* c"*02 

ture plane is at x = x^z which defines h' = cosa ' ^x * 

(y, z) coordinates where this ray strikes the aperture 
plane are: 

y ' ^02 <*c'^02^ cosa 'x ' ^02 ^*c"*02^ 

s. 


. / V cosa" z _ _ X rz2 

^ - ^02* <*o-’‘02> ^SSTx ' ^02 + <>‘c-*02' 


Then 


and the phase of the field upon reaching the aperture plane 


IS given as: 


2n 


^2 " ^ (R+RM+D) + Initial Phase. 
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Thus/ five par^uneter8 are computed for each ray at a point 
on the aperture pleme: the y and z coordinates, the y and 
z components of the electric field, and the phase of the 
field. 

B) Single Reflector System 

In this case each ray is traced from the feed to the 
reflector up to the aperture pleme in the same way as before. 
It is clear that in this case a smaller number of equations 
have to be written and the phase is given by 

2 IT 

I/) * — (R+d) + Initial phase. A more de- 

tailed discussion of the above operation is provided by 
Kauffman [ 2 ] . 

2.2 Calculation of Radiation Patterns 

In both cases, the tangent aperture field is given by: 


* (y + z E^„) e”^*^ for a single reflector 
AP ry rz 

where are the tangential components of the aperture 

electric field, or * (y E ^^2 * ®rz2^ * 

dual reflector. 

In order to evaluate the secondary radiation pattern 
at a particular point in space, we integrate numerically 
over the aperture. The integrals to be evaluated are; 


= // E cosd. sin(^ + z cose] 

. j rz 


Aperture 

Surface 


d d 

y z 


emd 
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B,= Fe sin0 + E COS0 

^ Aperture ^ •* 

Surface 

gjk sin0 siiKji + z cosoj^ ^ 

where the aperture surface is the area of the reflector 
aperture projected on the aperture plane. It is necessary 
to integrate only those points which result from reflections 
from the actual surface and not from its mathematical exten*> 
sion* This is achieved by interpolating a series of edge 
points on the boundary, using information from points which 
exist outside the aperture. All points then existing out~ 
side the reflector surface are disregarded. 

Before the integration takes place, all points on the 
aperture plane are quantized in their y-coordinate . All 
details on quantization and integration are fully provided 
by Kauffman [ 2 ] , Agrawal [ 3 ] , and Botula ("s] . 

2.3 Transition from the Old Algorithm to the New One 

The block diagram in Figure 2.9 shows the locations 
where changes, additions and modifications were applied to 
the old algorithm to obtain the new one. 

These general additions and changes, which will be ex- 
plained later in more detail, are the following: 

1. NPUT: Was enlarged to read in and print out 

data for both reflectors for a dual 
reflector antenna system. This feature 
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2 , SUBPNT: 


3. APRTUR: 


4. FIliDXC: 


5. CASSA: 


did not exist before. MPUT also calls 
an additional subroutine, named SUBPNT* 
Was added to determine the four extreme 
points on the subreflector, given the 
four extreme points on the main reflec-^ 
tor. 

Was extended for the following reasons: 

A) to incorporate hyperboloidal surfaces, 
as an addition to the previous list of 
surfaces . 

B) To ccxnpute, automatically, the loca- 
tion of the aperture pleme (x^) in 
terms of par^uneter8 pertinent to the 
antenna under consideration* This is 
accomplished by calling the 8\;U3routine 
FINDXC. 

FINDXC was added to provide APRTUR with 
an approximate value of x . x is eval- 
uated for both reflector systems, follow- 
ing different approximations depending on 
whether the antenna is a dual or a single 
reflector system. 

A new subroutine was inserted in APRTUR 
to account for all the tracing from the 
subreflector to the main reflector, up to 
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the aperture plane for the case of a 
dual reflector system. 

The rest of the program is unchanged. 



Fig. 2.9. Structure of new algorithm 
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3. STRUCTURE OF REFLECTR 


3.1 N»w Variables 

New variables were Introduced to account for the in- 
creased complexity of the program. Some old variables and 
conmon storage blocks were changed to give the new algorithm 
a general character. Since the new variables come an a 
follow-up of the old ones/ all common storage blocks and 
variables are introduced here. 

1} BLOCKG/YCBlJ ZCBL, BFMABL, HFMIBL (Aperture plane 
blockage information) . 

YCBL, ZCBL: y and z center coordinates of the aperture 
plane blockage ellipse. 

HFMABL, HFMIBL: 


2 ) 


Half-major and half-minor etxes of the aperture plane 
blockage ellipse. 


CASS/SR(3) , XO 
ER2(3), ER(3) 
XO, YO, ZO, 

X02, Y02, Z02 


V, Z 


RM 


, YO, ZO, Y, Z, RM, D, X02, Y02, Z02, 
(Only for Cassegrain antennas) . 

A point where a ray emanating from the 
feed intersects the subreflector. 

A point of intersection of the main 
reflector and the ray. 

The y and z coordinates of each ray on 
the aperture plane. 

Length of a ray from the subreflector 
to the main reflector. 
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3) 


4) 


D Distance of aperture plane from main 

reflector. 

SR (3) A vector s^^in the direction of a ray 

reflected by the subreflector. 

ER2(3) The three components of the electric 

field reflected by the main reflector. 

ER(3) The three components of the electric 

field reflected by the subreflector. 

COLOS/DELT, XC, ANGING, PM(3,4) , RS, XMX, ZMX, ZMN, 


YMX (Parameters used for determining x^.) 


DELT The 0" angle subtended by the subreflec- 

tor. (See Figure 2»2«) 

Angular increment. (See Botula for 
more details. 

Four extreme points on the main reflector. 
Distance from an extreme point on the sub- 
reflector to the origin. 

XMX, YMX, ZMX A point on the subreflector which is the 
closest point to the origin. 

ZMN The minimum Z coordinate of the subre- 

flector. 


ANGING 

PM(3,4) 

RS 


C0NTRL/N0PT(3) , NLIST, lOPT, ICASS, ILIST (100) 

NOPT(3) Three number specifying options regard- 

ing printer, plotter, and aperture plane, 
data output, respectively, i^ee [s] Section 
6 .) 
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The number of panels for which the algo- 
rithm will print complete illumination 
and quantizing data. 

A variable which is zero when the program 
is to run normally, and one when the 
single-panel option is in effect. 

A variable which is one if a Cassegrain 
antenna is to be analyzed, and zero for 
a single reflector antenna. 

The specifi<? panels for which the algo- 
rithm is to provide complete illumina- 
tion and quantizing data. (See Botula 
[sj , Section 4 . ) 

5) The common blocks; A) DIMENS, B) EXTENT, C) MATH and 
D) PATTRN, have remained the same as in j^sj. 

6) FEED/EIM91), ET(91), NP, N1 , XS, YS, ZS. 

(Peed antenna paraniotors) 

EP(91) , ET(91) 

Array containing the electric field strengths of 
the feed antenna in one-degree increments off-axis 
in the 0 - 90^^ and « 180^ planes, respectively. 
NP, NT 

The numbox' of increments of phi and theta used in 
the illumination pattern, respectively. 

XS, YS, ZS 

A point on each panel which is the closest point to 
the origin of the reference coordinate systtxm. 


MLIST 


lOPT 


ICASS 


ILIST(IOO) 



23 


PARAMS/AORORF , BELLP, CELLP, DIST, PSI, PLNPMT (3), 
PLNORM (3), FEED (3), ALPHA, BETA, GAMMA, XLAM, 

AOROR2, BELLP2, CELLP2, PSI2, DIST2, POINT (3), NORM 
(3), SURFCl, NPNL, NPOINT, SURFC2. (Antenna system 
parameters . ) 

In the following, the variables that appear first are 
defined on the subreflector, and those that appear second 
are defined on the main reflector. 

AORORF, AOROR2: The focal length of a paraboloidal re- 

flector, the focal length of a parabolic cylindrical reflec- 
tor, the radius of a spherical reflector, the semi-major 
axis of an ellipsoidal reflector along X, or half the trans- 
verse axis (x-direction) of a hyperboloidal reflector (Appen- 
dix B) , depending on which surface is intended to represent 
the reflector. 

BELLP, BELLP2: The semi -minor axes (along y and z, 

respectively) of an ellipsoidal reflector surface. Note 
that this does not define a completely arbitrary ellipsoid 
since the axes along y and z must be equal. For the case 
of a hyperboloidal reflector surface, this value represents 
the y semi-axis of the ellipse in the yz plane of the hyper- 
boloid. 

CELLP, CELLP2: Used only for a hyperboloid and stands 

for the z semi-axis of the ellipse in the yz plane of the 
hyperboloid . 
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DIST: A parameter used In translating the origin of 

the hyperbolic subreflector coordinate system so that it 
coincides with that of the main reflector. (See Appendix 
B.) 

PLNPNTO), POINT (3): The coordinates of a point on a 

planar reflector surface (x, y, z) . 

PLN0RM(3 ) , NORM (3) : The components of a unit normal 

vector to a planar reflector surface (x, y, z) . 

FEED (3): The reference coordinate system origin as 

expressed in the feed coordinate system (x, y, z) . 

XLAM: Wavelength of the feed antenna radiation. 

ALPHA, BETA, GAMMA; Rotation angles mentioned before 
(Figure 2.4.) 

SURFCl, SURFC2; Integer variables which determine the 
type of reflector surface. (This code is applied to the 
subreflector as well as the main reflector.) 

1) Surface is a plane. 

2) Surface is an ellipsoid. 

3) Surface is a sphere. 

4) Surface is a paraboloid. 

5) Surface is a parabolic cylinder. 

6) Surface is a hyperboloid. 

NPNL: Determines the number of panels the reflector 

is made of. The value of one means that a list of perimeter 
points and other surface parauneters for each panel must be 
supplied. In this case, the aperture boundary is approxi- 
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mated by a polygon. The value of zero means that the 
single-panel option is in effect and hence an ellipse 
is used to represent the boundary of that panel. 

NPOINT: The number of rays stored for processing in 

the P array at any given time. 

3.2 NPUT 

This is an input/output routine. If ICASS - 0, the 
program is to analyze a single reflector antenna system 
with two options: 

1) With lOPT = 1 for a single-pemel option. 

2) With lOPT = 0 for a multipanel option. 

In both cases, the four extreme points of the reflector 
surface are required. If ICASS « 1 a dual reflector anten- 
na is to be analyzed. For this case, the four extreme points 
of the main reflector are read in and used to find the four 
extreme points of the subreflector by calling subroutine 
SUBPNT. (SUBPNT explained later in this Section.) 

NPUT also reads other parameters concerning the feed. 
This is important since all pieces of information read here 
are used in conjunction with the FILL routine which is called 
later in the program. The connecting agent in this operation 
is the common storage block, named FEED. 

Previously, the four extreme points on the reflector 
were read into the P array only when NPNL was zero, and the 
variable x was also provided by the user. In this algorithm 

w 

the four extreme points are read regardless of the particular 
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value of NPNL. The reason for thi'i is that the above points 
are needed to compute the variable x later in the program. 
Furthermore, nevr printing statements were added to be useddd 
for dual reflector antennas. 

3.3 SUBPNT 



Fig. 3.1, Finding the four extreme points 

SUBPNT is called only fora dual reflector antenna . There 
is a "Do" loop which computes the distance (RR) from the ex- 
treme point on the main reflector to the reference point. 

RR = [^(PM(1,K))^+(PM(2,K) )^+(PM(3,K) )^] ^ 

vAiere PM(1,K) , PM(2,K) , and PM(3,K) are the coordinates of 
each extreme point on the main reflector. Then, the direc- 
tion cosines are found as: 

DIRl = PM(1,K)/RR (direction cosine in the x-direction) 
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DIR2 ■ PN(2,K)/RR (direction cosine in the y-direction) 

0IR3 ■ PN(3,K)/RR (direction cosine in the z-direction) 

The parametric equations of a line passing through the 
origin (reference point ) , and a point on the main reflector 
are given by: 

P(1,K) » PM(1,K) - RR«DIR1 

P(2,K) « PM(2,K) - RR*DIR2 

P(3,K) = PM(3,K) - RR»DIR3 

where P(1,K), P(2,K)/ P(3,K) is an extreme point on the 
reflector. To determine this point , the above parametric 
equations and the equation of the surface of the subreflec- 
tor are solved simultaneously. (See Appendix C for details.) 

This operation is repeated four times, i.e., once for 
each extreme point of the subreflector. 

3.4 APRTUR, APRIN, AND FILL 

APRTUR does all the ray tracing for the single reflec- 
tor antenna and it calls a new subroutine named CASSA for 
additional tracing in the dual reflector case. Figure 11 
shows the difference in approach between the old and new al- 
gorithms in determining the location of the aperture plane 
before integration for a multipanel, single reflector anten- 
na. 

This difference gives some increased accuracy in pre- 
dicting the radiation pattern of a multipanel, single reflec- 
tor emtenna. (See results. Section 5.) In the caue of a 
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single reflector, a short "Do" loop is used to find XHX, 
YMX, and ZMX, a point of the reflector which is the closest 
one to the origin. 



Then a rotation matrix A is computed from the rotation 
angles ALPHA, BETA, and GAMMA. The inverse of that matrix 
is also found. If the dual reflector option is in effect, 
the rotation matrix is calculated immediately skipping the 
above-mentioned "Do" loop. For single reflector antennas 
comprised of a number of panels, subroutine APRIN is called 
to provide data for each panel individually. 

Two important additions have been made in APRIN: 1) 

For each plane reflector a normal is computed automatically 
using the principle of the CROSS product. (See Appendix D.) 
2) Statements 20-28 make use of a "Do" loop to search for 
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(x / y f Z-) f A point on each panel, which is also the closest 
8 8 8 

point to the origin. It is an important point because it is 
used later, in APTRUR, to find the location of the aperture 
plane for each panel individually. (See Figure 3.2 for 

ge<xnetry) . For a complete discussion of APRIN, see 

From statements 50 to 65, APRTUR finds the angles sub- 
tended by the reflector or the reflector panel. Notice that 
in the dual reflector case, the angles subtended by the sub - 
reflector are the ones to be measured emd not those for the 
main reflector. All points, either the perimeter points for 
a pemel, or the four extreme points for a single panel option, 
are expressed as angles in the feed system. Then, a search 
for the maximum and minimxim 0' amd angles represented by 
the above-mentioned set of points is performed to determine 
the angles subtended by a panel or a subreflector. (See Ap- 
pendix B in [s].) 

ILLUMINATION ARRAY - Statements 65-95 generate the ap- 
propriate illumination array to insure a well-ordered illu- 
mination of the chosen reflector option. The previous method 
of illumination has been kept the same since it serves the 
purpose of the new algorithm in a rather convenient way. 

(See Section 2.3 in 

For the dual reflector case, the arngles subtended by 
the subreflector are the ones to be considered instead of 
those of the main reflector. The reason for this is the fact 
that an overillumination of the subreflector results in an 
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overilluxnlnation of the main reflector. Overlllvonination is 
desired so that the projected boundary of the main reflector 
on the aperture plane can be defined before integration is 
performed. The rays corresponding to the upper and lower 
limits of 6" miss the real subreflector. They get reflected 
by its mathematical extension, and as a result, they miss 
the main reflector too. 

If a Cassegrain antenna is to be studied, as soon as 
ANGINC is computed in APRTUR, subroutine FINOXC is called. 
(See Section 3.5.) This is the first time where FINDXC ap- 
pears in the program to provide APRTUR with the location of 
the aperture plane (x ) . APRTUR, with a "Do" loop in state- 

C 

ment 95, loads all illumination angles into the P array just 
after the angle pairs corresponding to the perimeter points. 
SUBROUTINE FILL is called to provide the angle pairs in the 
P array with the field strength and phase values. 

FILL - This routine is changed and adjusted to each 
antenna whose radiation pattern is to be computed. A detail- 
ed description of this subroutine and its various forms ap- 
pear in [ 2 ] , [^3j and [*5j . A new subroutine has been written 
for a vertical polarization case. (See Appendix E.) 

Furthermore, in APRTUR for single reflector antennas as 
the Zq) point is found, the location of a separate 

plane are determined. This part of the algorithm is not 
carried out for dual reflectors. The procedure for determin- 
ing XX and x^^ is as follows: 
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If the single panel option is in effect, then subrou- 
tine FINING is called. This is the second location in the 
program where FINDXC appears. (See Section 3.5.) If a 
multipanel option is in effect, then (Figure 12) is ex- 
pressed as: 

Rl » |(Xg + B12)^ + (Yg + B22)^ + (2g + B32)^‘* - 1.0 « R'-l.O 

where R* « jlxg + B12)^ + (y^ + B22)^ + (z^ + B32)^J** is the 
distance between (x , y , z ) and (B12, B22, B32) . 





Fig. 3.3. Computation of and xx 

Also, the angle subtended by the reflector is expressed 

as: 

®max * tan~^|- j ^ + 3 ^ 2 ) is a 

negative value and (z + B12) a positive one. Hence, to 

s 

obtain a positive angle, a negative sign is added. The 
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angle is augmented by 2.5 ANGIMG« i.e., 2.5 times an 

angular increment. The reason that * R' - 1.0 and 

9...^ “ 9.».„ + 2.5 ANGING are used Instead of R* and 9*.„» 
au 9 max max 

is to make sure that the pemel will be overilluminated. 

Thus is found as: 
c 

The distance between x and x for the first panel is com- 

c s 

puted as CONST ■lx - x | . This number beccxnes an impor- 
tant factor in locating the aperture plane for the rest of 
the panels. The idea is to put an aperture plane in front 
of every panel and with a distance equal to CONST away from 
it. This results in having an ordered arrangement of aper- 
ture planes in front of the reflector. So, the rest of the 
x^*s are given as: 

X = X + CONST where x„ is provided by APRIN, in ad- 
c s s 

vence. Once all x *s have been found, the location of a 

c 

general plane (xx) is determined, using FINOXC. (See Sec- 
tion 3.5.) Each panel is first projected onto its own indi- 
vidual aperture plane, and then phase-referenced to the gen- 
eral aperture plane. Thus, the general plane sums up all 
these projections that comprise the total projection of the 
2 mtenna on the aperture plane. This method of preparation 
of the aperture plane before integration yields better re- 
sults compared with the previous method. 

The difference in phase is written as: 
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DIFF - I X - X. I and tho I HASE - ~ (R+D+DIF) Initial Phase 

I C 8 I A 

where R ■ distance from the feed to reflector* 

D ■ a distance from the reflector to the Individual 
aperture plane. 

DIF « distance from the individual aperture plane to 
the general one. 

If the dual reflector antenna option is in effect, sub- 
routine CASSA is called by APRTUR to continue the ray trac- 
ing operation over the region lying between the subreflector 
and the main reflector* (See Section 3*6.) 

3.5 FINDXC 


This subroutine is called, as mentioned before, at two 
different locations in APRTUR. 



Fig. 3.4. Location of an aperture plane at for a 

dual reflector 
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In the dual reflector antenna case, FINOXC is called 
immediately after AMGING is confuted. In this case, x_ is 
evaluated directly from the geometry of the two reflectors. 
From Figure 3.4, a point with the largest z coordinate on 
the main reflector is determined and its distance (R') from 
the reference system is computed. Then, new parameter RSM 
is c<miputed as: 

RSM « [(PM(1,M))^+(PM{2,M))^•^(PM(3,M))^J**«1.0 • R' - 1.0 
Where R'« ^ (PM(1,M) )^-MPM(2,M) )^+(PM(3,M) ) 

Also, the angle subtended by the main reflector is 

niAX 

expressed : 

— 1 PM (3 M ) 

®max “ ^’“ PMtS ^M) ^ negative sign is 

provided here to obtain a positive angle, since 

PM(3,M) is positive and PM(1)M is negative. In the refer- 
ence system another angle, called 6 augmented is estimated as 

6aug « 6max >»■ 3.0 * ANGINGdn radians) and x. is then 

c 

calculated using the expression. 

x^ * -RSM cos (6aug.) 

The fact that RSM is used instead of K* and 

instead of is to insure overillumination and to make 

inAx 

sure that this subroutine works for all sub and main reflec- 
tor combinations, no matter what their geometrical relation- 
ships are. This subroutine could, if necessary, be changed 
to deal with each sub and main reflector c< ' binations s*?p- 
arately . 
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The second call of FINOXC by APRTUR is concerned with 

finding the location (xx) of the general aperture plane for 

the multipemel option reflector, or x_ for the single pemel 

c 

option. This task is accomplished as follows: (See Figure 

3.5.) 

✓ 

✓ 



Fig. 3.5. Location of xx or x^ for multipanel or single 

panel antenna, respectively 


First, find the distance R' between the point 
Zmx) , and the feed, i.e.. 


R' = 


(x^^ + B12) 
mx 


* ® 22 ) 


+ (z, 


mx 


^*mx' ^mx' 

+ B32)^l^ 


where y^, z^) is the point on the reflector which 

is the closest to the origin of the reference system. It 
should be noted that this point is computed at the beginning 
of the APRTUR routine. Second, 


3C 


e 


max 


tan 


-1 


^ ® 32 ^ 


^ I gives the maximum angle 


^nx ' “12' 

subtended by the reflector. This angle is increased by a 
2.5 AMGING to give 6aug * Omax -f 2.5 ANGING (in radians) 
and third, to find x , R' is reduced by 1.5 to yield 


RSM 
and hence 




+ (y„v + 

mx '*mx mx i2 


- 1 .! 


x^ or XX = - RSM cos(Oaug) + B(l,2) for a single 
panel or a multipanel antenna, respectively. 

It is noted here that the distance R' is reduced by 

1.5 instead of 1.0 (as was done in the case of individual 
panels) to insure that xx will be less than x , in the multi** 
panel case. The whole arrangement of separate aperture 
places and a general one is shown in Figure 11, Part B. 

It can be seen that xx has to be behind all individual 

aperture planes. If the multipanel option is not in effect, 

XX becomes x^. 

c 

3 . 6 CASSA 


This subroutine accomplishes all the ray tracing from 
the subreflector to the main reflector up to the aperture 
plane. It starts with finding the direction cosines of a 
vector along the ray reflected by the subreflector. Para- 
metric equations of a line are expressed as: 

Xq2 = ^0 BM*DC(1) 

yo2 = yg bm*dc(2) 

2q2 = *0 BM*DC(3) 
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where (Xq 2 # Yq 2 * ^ point on the main reflector. 


(Xq« Yq* Zq) is a point on subreflector, RM distance be- 
tween these two points and DC(1) DC (2) DC (3) are the direc- 
tion cosines with respect to x, y and z axes, respectively. 
The solution of simultaneous equations consisting of the 
above parametric equations and the equation of the reflector 
surface yield the point Xq 2 , Yq2* *02* Although this sub- 
routine has been written to deal with six analytical sur- 
faces, it could be extended to incorporate any other number 
of types of surfaces, if desired. Surfaces expressed numeri- 
cally could also be added to this algorithm, especially for 
the dual reflector antenna option, where shaping of one or 
both of the reflectors is now widely used in their actual 
design. 

Once the point Xq 2 » yQ2' ^02 evaluated, the normal 
(NHAT2(1), NHAT2(2), NHAT2(3)) on the surface at that point 


is computed as follows: 

Let the surface be represented as g(x, y, z) 

Vg(xQ2» Yo2' *02^ 


= C. 


A 

Then n 


02 




A detailed explanation of computing normals and intersec- 
tions of rays with surfaces is not given in this thesis, 
since a complete discussion can be found in all references 
from ^ij to in their description of subroutine APRTUR. 

The only difference lies in the fact that the parameters used 
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in CASSA are pertinent to the surface of the main reflector 
and not the suhref lector. 

The normal on the main reflector is used to apply 
Snell's law of reflection to find a vector in the direction 
of the reflected ray (SR2(1), SR2(2), SR2(3)). This part of 
the algorithm is described in Section 2.1. A point, (y, z) 
on the aperture plane is then computed, and passed over to 
APRTUR where it is stored, t** he retrieved later by QUANTZ. 

The principles of geometrical optics are used to deter- 
mine the electric field during these two phases of ray trac- 
ing. All equations in this part of the algorithm are men- 
tioned in Section 2.1. In general, all operations taking 
place in CASSA arc depicted in Figures 2.6 emd 2.8. 

3.7 Main Procedure and the Utility Routines 

The main procedure and all the rest of the utility 
subroutines were kept the same as before with a minor change 
in their storage blocks. A complete development of these 
subroutines and the main procedure is provided by Botula in 

W- 
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4. EXAMPLES AND TEST CASES 

4.1 Introduction 

Two test cases on the Cassegrain antennas are provided 
here to demonstrate the use of the program and support the 
validity of the algorithm. These cases are the following: 

FIRST / a classical Cassegrain antenna which was used 
to check the algorithm in the case of uniform illumination, 
but with no blockage. 

SECOND , a dual offset reflector antenna, used to check 
the results obtained by this algorithm against clculated data 
obtained from two other algorithms. 

4.2 Exeunple and First Test Case 

The classical Cassegrain antenna, shown in Figure 4.1 
employs a hyperboloid for the subreflector and a paraboloid 
for the main reflector. One of the two foci of the hyper- 
boloid is the real focal point of the system, and is located 
at the origin of the feed coordinate system; the other is a 
virtual focal point which is located at the focus of the 
paraboloid which coincides with the origin of the reference 
system. As a result, all rays originating from .the real 
focus and reflected from both surfaces travel equal dist^mces 
to a plcuie in front of the antenna. (See Figure 4.1.) 
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PARABOLOID HYPERBOLOID 



Table 4.1 gives a number of parameters that define 
completely the geometry of the antenna system. All para- 
meters required by NPUT will now be evaluated from this 
Table . 

TABLE 4.1. 

Main reflector focal length (F ) 

m 

Main reflector illumination angle (Oj) 

Eccentricity of subreflector (e) 

Distance between two foci (F ) 

Wavelength (X) 


100.0 in. 
60.785 
1.50177 
91.0 in. 
4.734 in. 
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FEED PARAMETERS 

Since the origin of the feed coordinate system is 
located at the real focus of the hyperboloid and at dis- 
tance X = -91.0 from the origin of the reference system, 
the feed parameters can be given as: 

1) Feed (1) = -91.0 in.. Feed (2) =» 0.0, Feed (3) =0.0 

2) ALPHA = 0.0, BETA = 0.0, GAMMA = -180.0 
MAIN REFLECTOR PARAMETERS 

SUBFC2 is set equal to 4, since a paraboloidal reflec- 
tor is to be used as a main reflector. 

F = 100.0, as was given in the Table, 
in 

The four extreme points of that reflector to be read 
in are; 

Upper point 

r. - ^ ^m _ ^ ^m , 2 .( 100 . 0 ) = . 04 4 

“ 1 + COS e * 1 4- cos 6 ^ 1 + COS ( 6 6. 1 6 ^ ) 

max 

X = r cos = -r cos(60.785°) = -65.599 

max 

y = 0.0 

z = r sin = r sin(60.785°) = 117.304 

max 

Lower point 

r- - ^^m _ 2.100.0 - A A 

” 1 + cos ” 1 + cos(-60.785^J ” 
mxn 

X = -r cos “= -65.599 

mm 

y = 0.0 

z = r sin = r sin (-60.785°) = -117.304 
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These two points correspond to the Q' extrema in the feed 
system. Also, the two points representing the y - extrema 
are almost exactly the extrema as well. The z coordi- 
nates of these points are identical. 

_ _ *min *max _ 117.304 - 117.304 

a B 2 2 0.0 

The reflector, as seen from the geometry of the antenna 
system, is 234.608 inches wide and symmetric with respect 
to the xz pleme, hence y = = ±117.304 in. 

Finally, the paraboloid equation provides the x coordinates 

X . - “.599 

m 

Thus, the four aperture points become: 

Upper point: (-65.599,0.0,117.304 = 

PM(1,1),PM(2,1) ,PM(3.D) 

Lower point: (-65.699,0.0,-117.304 = 

PM(1,2),PM(2,2),PM(3,2)) 

Leftmost point: (-65.599,-117.304,0.0) = 

PM(1,3) ,PM(2,3) ,PM(3,3)) 

Rightmost point: (-65.599,117.304,0.0) = 

PM(1,4) ,PM(2,4) ,PM(3,4)) 

It should be noted here that the diameter of the main reflec 
tor can also be found from the relationship given in Appen- 
dix A as follows: 


tan 


1 



+ D 


m 


4F^ tan 

m 


60.785° 

5 


= 234.608 
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SUSREFLECTOR PARAMETERS 

SURFCl is set equal to 6, since a hyperboloidal surface 
is to be used for a »ubref lector. NPNL takes the value of 
zero, since neither the subreflector nor the main reflector 
is composed of panels. 

The parameters a (semi-transverse axis along x > 

AORORP) , 

b (semi-axis along the y direction = 
BELLP) , and 

c (semi -axis along z direction « CELLP) 
are computed as follows: (See Appendix B for details.) 

^ “ "T * rXSOI?? “ 30-2976 

c = b = a ^e^-1 = 30.2976 | (1.50177)^-1 » 33.95 

F 

Also, DIST ^ ~ = 45.0 which is a parameter used in 

translating the origin of the subreflector coordinate system 
so that it coincides with that of the main reflector. (See 
Appendix B for details.) 

There is no need to read in x_, since this value is 

c 

computed in the program as a function of the emtenna system 
parameters. In this case, the FILL routine was not used, 
and no data for the E and H plane patterns of the feed were 
used in the input file. 
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4.3 General Input File 

For format information. 
Appendix F. 

A) Dual Reflector Cases 
Cards 
1-4 

5 

6 

7 

8 

9 

10,11,12,13 

14 

15-N 

N+1 

N+2 

N+3 


refer to the program listing. 


Information 
Title Cards 

Feed (1-3) , ALPHA, BETA, 
GAMMA, XLAM 

SURFC2, AOROR2, BELLP2, 
CELLP2, DIST2, PSI2 

P0INT(l-3), NORM(l-3) 

SURFCl, NPNL, AORORF, 

BELLP, CELLP, DIST, PSI 

PLNPNT(l-3) , PLNORM(l-3) 

Four extreme points 

edge of '^rhe'^main reflector. 
One point goes on each card. 

YCBL, ZCBL, HFMABL, HMIBL 
(Blockage of main reflector 
by subreflector) 

Any data required by the 
FILL routine 

NOPT, NLIST 

MAJOR, AMAJOR, MINOR, 
AMlNOR(l-3) (Pattern request 
cards) 

DONE typed in the first 
four columns of the card 
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B) Single Reflector Cases 
Cards 
1-4 

5 

6 
7 

8,9,10,11 

12 

12-N 


Information 
Title Cards 


Feed (1-3) , ALPHA, BETA, 
GAMMA, XLAM 


SURFCl, NPNL, AORORF, 
BELLP, CELLP, DIST, PSI 


PLNPNT(l-3) , PLNORM(l-3) 

Four extreme points 

(x, y, z) on the reflector. 

One point goes on each 

card (single panel option 

only) 

YCBL, ZCBL, HFMABL, HMIBL 
(Blockage of reflector fy 
ieed) 

Any data required by FILL 
ROUTINE 


N+1 NOPT, NLIST(if NOPT speci- 

fies that only certain 
panels are to be printed or 
plotted , cards containing 
the list of these panels 
follow this card) 


N+2 MAJOR, AMAJOR, MINOR, 

AMINOR(l-3) 

N+3 DONE is typed in the first 

four columns of the card 


These cards are followed by the panel data 
of the panel data is as follows: 


The organization 
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Information 

NPERIM, SURFCl, NPTPPL 

(x, y/ z) perimeter point 
(one perimeter point per 
card) 

AORORF, or AORORF and 
BELLP, or AORORF and PSI, 
or PLNPNTO), or AORORF, 
BELLP, CELLP, and DIST, 
depending on which para- 
meters are needed to des- 
cribe the surface speci- 
fied by SURFCl. 

All cards carrying information for individual panels appear 

in the main input file after the DONE card. 

4 .4 Development of a Uniformly Illuminated, Classical 
Cassegrain Antenna 

All parameters needed for this case were computed in 
Section 4.2. None of the available FILL subroutines was 
used emd the H and E plane patterns (for the feed) were not 
read in as data in this particular case. The reason for 
that was to insure uniform illumination over the main reflec- 
tor. The procedure adopted to achieve this task was as 
follows: 

1. FILL is not called in APERTUR. 

2. All lines in APERTUR related to the 
amplitude and phase of the E field were 
moved to subroutine CASSA. 

3. In subroutine CASSA the following modifica- 
tions took place: 

PI = R‘RM and P2 = 0.0 


Cards 

1 

2-(NPERIMtl) 

NPERIMt2 
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®Ti “ IT equal to RM) and “ *R" * 0.0 

where and Ep^ are the e and ^ electric field con^nenta 
of the incident (on the aubref lector) ray, reapectively. 


From end applying Snell 'a law to raya reflected by the 

two aurfacea, and E ^2 ^ere evaluated, where E^ ia the 
electric field vector along a ray reflected by the aubre- 
fleeter, and £^^2 electric field vector along a ray 

incident on the main reflector. 

It ia obvioua that in the far field, £^^2 * ^ 


E^i/R 

Also, 


Pl/R , PI 
RM R*RM 


R«RM 

R*RM 


1.0 


Which means that th^^ £ field waa kept conatant at the value 
of one along every ray. Thua, the conatant amplitude re- 
quirement for uniform illumination was met. 

4. The constant phase requirement waa also 
satisfied by the above arrangement, since 
the phase was set equal to: 

PHASE « ^(R+RM+D) . 

Notice that R+RM+D is always constant for 
a focused Cassegrain antenna. (See Appendix 
A.) 

Table 4.2 shows the input file for Case A. The first 
four cards contain title information which is also reproduced 
at the printout. Information about the feed coordinate 
system (FEED, ALPHA, BETA, GAMMA, and XLAM) appear on Card 5. 
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Cards 6 and 7 contain Infomation about the surface of the 
main reflector. Card 6 is for SURFC2, A0R0R2, BELLP2, 
CELLP2, DIST2, PSI2 and Card 7 Is for POINT, NORM. For 
this main reflector, SURPC2 » 4 and A0R0R2 s 100* C. None 
of the other parameters is required for this surface, so 
all are given the value of zero. Cards 8 and 9 contain 
infomation for the subreflector surface. Card 8 is for 
SURFCl, NPNL, AORORF, BELLP, CELLP, DIST, PSI and Card 9 
for PLNPNT and PLNORM. For that type of subreflector sur- 
face SURFCl « 6, NPNL - 0.0, AORORF « 30.2976, BELLP « 

33.95, CELLP « 33.95, and DIST « 45.500. The rest of the 
other parameters are given the value zero, since none of 
them is required for this surface. Cards 10, 11, 12 and 13 
contain the four extreme points (x02, Y02, Z02) of the main 
reflector. Card 14 carries the required blockage informa- 
tion, i.e., YCBL, ZCBL, HFMABL, and HFMIBL. In this case, 
aperture blockage is not considered and so all the above 
parameters are set equal to zero. Since the FILI< routine 
is not used in this case and no data for the feed radiation 
patterns are needed. Card 15 is used to deteimine the output 
option code. Here the computer is instructed to print and 
plot information about the two surfaces, as follows: 

NOPT(l) 2 (print all results) 

NOPT(2) s 2 (plot apertxure after quantizing) 

N0PT(3) - 1 (print aperture array onto a disc 

file at the end of QUANTZ) . 
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NLIST is equal to zero since the antenna in question is not 
divided into panels. Cards 16 and 17 are the radiation pat- 
tern requests. One pattern is required in » 0^ plane for 
6 from 85.0° to 95.0° by increments of 0.5°, and another one 
in the 0 = 90. 0° plane for = -4.0° to 4.0° by 0.5°. The 
next and last card (No. 18) has DONE typed in the first four 
coltunns, which signifies the end of the pattern requests and 
the end of the input file. The result of this check case are 
shown in Appendix G. 

Figure 4.2 shows a comparison of the results obtained by 
this algorithm with those results reported by Silver for a 
uniformly illuminated circular aperture [oj . 


so 


UNIFORMLY ILLUMINATED 
CLASSICAL CASSEGRAIN ANTENNA 



Fig. 4.2. Classical Cassegrain antenna radiation pattern 
(Due to the syinmetry, only one-half of 
the pattern is shown) 


TABLE 4«2; 

CASE A INPUT FILE 

1 CASSEGRAIN ANTENNA EXAMPLE 

2 A PARABOLOID-HYPERBOLOID COMBINATION 

3 FEBRUARY 13, 1981, NCSU PGMR: CHRISTOS 

FCLTY;RD-DAN 
PRT: HILLSBORO 

4 (A BLANK CARO) 


5 

-91.U05 

0.0 

0.0 

o 

» 

o 

O 

« 

o 

-180.0 

6 

4 

100.0 

0.0 

0.0 

0.0 

0.0 

7 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

8 

6 

030.2976 

33.95 

33.95 

45.50 

0.0 

9 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

10 

-65.5997 

-117.304 

0.0 




11 

-65.5997 

117.304 

0.0 




12 

-65.5997 

0.0 

-117.304 




13 

-65.5997 

0.0 

117.304 




14 

0.0 

0.0 

0.0 

0.0 



15 

221 






16 

PHI 

0.0 

THETA 

85.0 

95.0 

0.5 

17 

THETA 

90.0 

PHI 

-4.0 

4.0 

0.5 

18 

DONE 







51 


4.734 


Here, the algorithm is tested with calculated data 
reported by TICRA A/S [s], and C. C. Chen [s] . The reason 
for choosing an offset case as a second test case is the 
fact that offset geometry does not have the symmetry of 
the first test case, which can sometimes mask errors. 
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TABLE 4.3 


e 

Pc 

X 

Fro 

Offset angle (0°) 
Aperture diaroeter 


2.47 

33.07 in 
0.98425 in 
69.685 in 
37.60 

(Dm) 64.8 X 
16.40 


Tilted angle of 
feed axis (ai) 

-11 db taper was used. 
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Using the relationships between the hyperboloid and 
paraboloid from Appendix B, and using the given data in 
Table 4.3, one can estimate AORORF, BELLP, CELL? and DIST. 
Furthermore, in this case, ALPHA « 0.0, BETA « 0.0 and 
GAMMA s -163.6, since the axis of the feed makes an angle 
ioLi) of 14.6° with the x axis of the reference system, as 
shown in Figure 4.3. Feed (1), Feed (2), and Feed (3), as 
well as the four extreme points of the main reflector are 
easily calculated. The input file is shown in Table 4.4. 

In this case, the input file is arranged in the same way 
as before up to the fourteenth card. Cards 15 to 52 contain 
information cd>out the feed radiation pattern. Card 53 con- 
tains NOPT, NLIST, and Cards 54 and 55 are used for the 
pattern requests. Finally, DONE is typed on Card 56. The 
secondcury radiation pattern is shown in Figure 4.4, and com- 
pared with data obtained from the other two algorithms. 
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TABLE 4.4 
CASE B INPUT FILE 

1 OFFSET CASSEGRAIN ANTENNA EXAMPLE 

2 A PARABOLOID - HYPERBOLOID COMBINATION 

3 FEBRUARY 19, 1981 NCSU PGMR-CHRISTOS FCLTY:RD-DAN 

4 TICRA AP/S 


5 

-31.725344 

0.0 9.337276 0.0 

0.0 

-163.6 .98425 

6 

4 

69.685 0.0 

0.0 

0.0 

0.0 

7 

0.0 

0.0 0.0 

0.0 

0.0 

0.0 

8 

6 

6.694507 15.119643 15.119643 16.535433 

o 

o 

9 

0.0 

0.0 0.0 

0.0 

0.0 

0.0 

10 

-57.18974 • 

-31.88699 49.66035 




11 

-57.18974 

31.88699 49.66035 




12 

-45.82776 

0.0 81.54734 




13 

-68.55171 

0.0 17.77336 




14 

0.0 

0.0 0.0 

0.0 



15 

1.000000 

.99053 

.96266 

.91793 

.85878 

16 

.78830 

.70997 

.62737 

.54392 

.46269 

17 

.38617 

.31623 

.25407 

.20029 

.15491 

18 

.U756 

.08755 

.06394 

.04563 

.03223 

19-32.00000 

.00000 

.00000 

.00000 

.00000 

33 

.00000 





34 

1.0000 

.99053 

.96266 

.91793 

.85878 

35 

.78830 

.70997 

.62737 

.54392 

.46269 

36 

.38617 

.31623 

.25407 

.20029 

.15491 

37 

J.1756 

.08755 

.06394 

.04563 

.03223 
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Table 4.4 (continued) 


38- 

51 .00000 

.00000 

.00000 

.00000 .00000 

52 

.00000 




53 

221 




54 

PHI 0.0 

THETA 

87.0 

93.0 0.25 

55 

THETA 90.0 

PHI 

-3.0 

3.0 0.25 

56 

DONE 





CALCULATED 



Offset Cassegrain E-plane pattern 
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5. A SINGLE REFLECTOR ANTENNA EXAMPLE 
(A SEGMENTED SPHERICAL REFLECTOR) 

5.1 Description of the Problem 

A single reflector composed of 54 panels was construc- 
ted and tested by NASA at the Langley Research Center. Its 
measured radiation patterns were compared twice: First, 

with calculated results obtained by using the old version 
[sj; and second, with calculated results obtained via the 
modified version incorporated in the new algorithm. A com- 
plete description of the antenna and its parameters is pro- 
vided by Botula in [5]. Here, the input file and the 
results only are given. 

5.2 Results and Comments 

Figures 5.1 and 5.2 depict the projections on all panels 
on the aperture plane. The result obtained by the old ver- 
sion is shown in Figure 5.1, whereas the result from the 
revised algorithm is shown in Figure 5.2. 

Figures 5. 3-5. 6, inclusively, show the secondary radi- 
ation pattern for both versions. The reason for this 
discrepancy in the above results lies in the amount of 
overlapping between the projected panels on the aperture 
plane. The more the overlapping, the less accurate results 
are obtained compared to measured data. 

The reason for this overlapping is due to the fact that 
the rays reflected by the perimeter points of each panel tend 
to diverge on their way to the aperture plane. To reduce 
their divergence, the aperture plcuie is brought closer to 
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each panel so that the rays travel over shorter distances 
before they strike the aperture plane. Once this occurs, 
the projected panel is then phase referenced to the general 
aperture plane. 

This procedure, which is summarized in Figure 3.2, 
yields less overlapping and better results than the old 
version. 

5.3 Input File 


TABLE 5.1 

INPUT FILE FOR A SINGLE REFLECTOR ANTENNA 

1 Faceted Spherical Reflector Test Case (LSST) 

2 Surface composed of 54 panels, three perimeter points 

per panel, 

3 no blockage. Feed phase center 0.5 lamda inside horn 

aperture. Emplane only 

4 (Blank Card) 


5 

9.441 

0.00 

8.026 

0.0 

0.0 

-40.0 0.3335 

6 

3 5424.0 

0.0 

0.0 

o 

• 

o 


7 

o 

• 

o 

0.0 

o 

• 

o 

0.0 

o 

• 

o 

0.0 

8 

-19.6617 

O 

• 

o 

13.7628 




9 

-20.7015 

5.0252 

11.0542 




.0 

-22.2326 

-5.0581 

7.4916 




.1 

-23.8206 

0.0 

2.9285 




.2 

0.0 

0.0 

0.0 

O 

o 




13-50 Illumination data for FILL routine 
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3 
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Table 5.1 (continued) 


52 

1 

15 51 



S3 

THETA 

91.0 

PHI 

0.0 15.0 0.5 

54 

PHI 

o 

• 

o 

THETA 81.0 100.0 0.5 

55 

DONE 




56 

3 

1 150 



57 

XI 

Y1 

Z1 


58 

X2 

Y2 

Z2 

Three points for the first panel 
on reflector 

59 

X3 

Y3 

Z3 


60 

3 

1 150 



61 

XI 

Y1 

Z1 


62 

X2 

Y2 

Z2 

Three points for the second pemel 
on reflector 

63 

X3 

Y3 

Z3 



64-273 This process is repeated for all 54 p^ulels. 


2-flXIS 
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MAP OF PANEL PROJECTIONS 

15.0 -1 
12.5 - 

10.0 - 

7.50 - 

m 

5.00 - 

m 

2.50 - 

m 
m 

-7.50 -5.00 -2.50 0. 2.50 5.00 7.50 

Y-flXIS 



Fig. 5.1. Old algorithm 


2-ftXIS 
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MAP OF PANEL PROJECTIONS 


15.0 ^ 

m 

12.5 - 

10. 0 - 

7.50 - 
5.00 - 

2.50 4-r-r 
-7.50 


^ I • » • “• j » f "f I I 'l*T 

-5.00 -2.50 0. 

Y-nxis 


• I » I I I I I I I 'f I 

2.50 5.00 7.50 



Fig* 5*2* New algorithm 



Fig. 5.3. Sphere Emplane pattern (old algorithm) 
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6. CONCLUSIONS 

An algorithm capable of computing radiation patterns of 
single reflector antennas has been modified and extended to 
analyze dual reflector antennas. A new technique for deter- 
mining the aperture plane for multipanel single reflector 
antennas has been incorporated into the new program. The 
location of any aperture plane and the normals on each plane 
panel are computed automatically. Furthermore, equations 
for hyperbolic surfaces have been added. 

The capability of expressing any non-analytic surface 
numerically will render the present algorithm very versatile. 
This fact will make the analysis of dual reflector antennas 
with shaped surfaces possible. 

Presently, the algorithm requires that the feed center 
coincide with the real focus of the hyperboloid for a Casse- 
grain antenna, but modifications could be inserted to deal 
with any off-focus applications. 

The results for the dual reflector antennas obtained by 
this algorithm show good agreement with those obtained by 
other algorithms. It is believed that a direct comparison 
with measured patterns will give a better estimate of the 
accuracy of the present algorithm. 
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APPENDICES 
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8.1. APPENDIX A 
CASSEGRAIN ANTENNA GEOMETRY 



Fig. 8.1. Geometry of classical Cassegrain antenna 

The classical Cassegrain geometry shown above employs 
a parabolic contour for the main reflector and a hyperbolic 
contour for the subreflector. One of the foci of the hyper- 
boloid is the real focal point of the system emd is located 
at the origin of the feed coordinate system; the other is a 
virtual focal point which is located at the focus of the 
paraboloid. As a result, all parts of a wave emanating from 
the real focal point and then reflected from both reflector 
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surfaces, travel equal distances to a plane in front of 
the antenna. 

Four fixed parameters are adequate to completely des- 
cribe a Cassegrain system, two for each reflector. In 
Figure 8 . 1 , seven parameters are shown. If four are known, 
the other three can be derived frcxn the mathematical rela- 
tionships between the two reflector surfaces. For the main 
reflector. 


tan j 

for the subreflector: 


A = 1 2m 
®2 I Fm 


and 


2 F 


tcin6 1 tan 6 2 


2U'id 


± — .. ^ Z rm — 






where: F^ - distance between two foci, 

fj^, ^2 ~ lengths of hyperboloid. 

Dm = diameter of main reflector, 

D = diameter of subreflector, 
s 

Fm = focal length of paraboloid 

02 = one-half of the angle subtended by the main 
reflector 

01 = one-half of the angle subtended by the sub- 
reflector. 

For example, if Dm, Fm, F^ and 0 i are determined by consider- 
ations of antenna performance and space limitations, then 
02, Dg, and £2 can be derived. 
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Note 6# which determines the beamwidth required of the 
feed radiation pattern, may be determined independe.-tly of 
the ratio Fm/Dm which specified the sape of the main reflec- 
tor. 


The surface of the main reflector is given by: 

2 2 

y + z = 4 Pm (Pm + x ) , and the surface of the subreflector 
is expressed as: 


(X 


+ DIST) 
~2 


-h- 


b- ? 


= 1 


where DIST = 


_ c _ 


= a+ 


-Xq I (S ee Figure A. 2) is the distance 
used to translate the origin of the hyperbola coordinate 
system so that it coincides with the origin of the referenced 



a = half the transverse axis (along x-axis) 

b = semi-cixis along the y direction in the ellipse 
lying in the yz plane. 

c = serai-axis along the z direction in the ellipse 
lying in the yz plane. 


If e (eccentricity) of the hyperboloid is known, the 


following equations can be used: 


e 


sin% 

sin^ 


(62+61) 

TTa^STT 


a 



b 


a 



- 1 , 



E + 1 


M 


where M is the magnification factor of the hyperboloid 
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8.2. APPENDIX B 
ADDITION OF HYPERBOLOID 

The eqfiation of the hyperboloid, depicted In Figure 
8.3, In the cartesian system Is given as: 

x2 y2 z2 

-V > -y > -V « 1, where a ■ half the transverse axis 

iT c^ along x 

b » semi-axis of the ellipse In 
the yz plane 

c » seml-eocls of the ellipse in 
the yz plane 


VIRTUAL FOCAL 

REAL FOCAL SUBREFLECTOR poiNT 



In this equation, the hyperboloid Is expressed In th3 x , y 

s s 

and Zg coordinate system. To express the same surface in 
the x, y and z system, a translation has to take place along 
the X axis, so that the origins of the two systems 0 and 0 

S 

coincide. It is clear that y. * y and z_ » z, emd hence no 

s s 
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charge la needed to he made in the y and z directions. 

If DIST la the distance between O. and 0, then x can 
be expressed as x > x • DIST, or x ■ x OIST, and hence 

B 8 

the hyperboloid equation in the x, y, z system becomes: 

l i L . ii . where DIST « ^ (8.1) 

a^ b' c^ ^ 

and Fc “ distance between 
the two foci of 
the hyperboloid. 

The parametric equations for a ray are: 


X 

y 

z « RB 


*®11 ■ ®12 


*®21 “ ®22 


( 8 . 2 ) 


31 * ®32 


Substitute Equation (8.1) back into the equation of the 
hyperboloid to obtain: 

(RBjj - . DIST)2 . (RBj^ - Bjj) ^ . (BBjj - Bjj)^ . j . Q 


(8.3) 


or 


"" . ®12 . DIST^ BjjBjj 2R B^i DIST 2BjjDIST 

— T2 — 

AAA A A A 


R^b|j B^j 2 R BjjBj 2 B^B^j b ^2 . 2B31B3J 

— 7 ? ::5 :t ~::5 — ^ “ 

b b b c c c 


Equation (8*3) is of the form 
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AR + B R -f C - 0 


(8.4) 


where 


B 


11 


B 


21 


®n 

■? 


B 


®11 ®12 


Bj^j^ DIST 


-2 

,(DIST)2 

— J + — 5“' 5 J 

a a a^ b 


®21 ®22 


®31 ®32 




(8.5) 


( 8 . 6 ) 


(8.7) 


Equations (8v5), (8.6), and (8.7) are evaluated by the pro- 
gram and (8.4) is solved to find the intersection point of 
the ray with the surface. 

Now, to find the inside normal of the surface, the 
gradient of Equation (8.1) is taken as: 

7g(x,y,z) » n(x,y,2) (8.8) 


where 


g(x,y,z) 


(x DIST)^ z^ 


. C-ij-l (8.9) 
b c 


it follows that: 

* 3x y 3y * 3z 


- 2(xt.DIST), 2y,- 2z 

b^ c^ 


or 



2(x > DIST) 



i£ 


ay 



is 


az 



( 8 . 11 ) 
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Normalization results in obtaining the unit vector n as: 

n = ^ jx,y,z) _ X 2(x + DIST)/a^ - y (2/b^) -z (2z/a^) 

I tu * DIST)^ , 4y_ , 4z^ P (8.12) 
^ b* *? / 

The factor 2 Ccuxcels out from both numerator and denominator. 
Let the dencmiinator be expressed as: 


DEN = 


(x + DIST) 

i 


2 2 
b c 


(.8.13) 


then 


n = 


(x t DIST) /a 
DEN 


n = 


_ -y/b" 


DEN 


n = 


_ -z/c 
DEN 
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8.3. APPENDIX C 
SUBROUTINE SUBPNT 



Fig. 8.4. Determination of subreflector four 
outermost perimeter points 

In this subroutine, the four extreme points of the main 
reflec<*or are used to find the four extreme points on the 
subreflector. This task is accomplished as follows: 

Take a given extreme point on the main reflector and 
write the parametric equations of the line (RR) connecting 
that point to the origin of the reference system (0) . 

Express the direction cosines as: 


DIRl = COSA = PM(1,K)/RR 

(8.14) 

DIR2 = COSB = PM(2,K)/RR 

(8.15) 

DIR3 = COSC = PM(3,K)/RR 

(8.16) 
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Hence, the parametric equation of that line is given 


Xq = P(1,K) * PM(1,K) - RR'DIRI 
Yq * P(2,K) = PM(2,K) - RR*DIR2 
= P(3,K) = PM(3,K) - RR‘DIR3 


(8.19 


where (P(1,K), C*(2,K) and P(3,K)) is a point on the sub- 
reflector which is to be found. 

Now, substitute Equations (8.17), (8.18), and (8.19) in 
the equation for the surface of the hyperboloid, that is in 

(X tDIST)^ - zi - = 1 (8.20) 

a^ b^ c 

to obtain: 

[PM(1,K) - RR-DIRl + DIST] ^ _ [PM(2,K) - RR-DIR2] ^ 

_ 

a b 

rPM(3,K) - RR>DIR3]^ _ , ,, , 


(PM(1,K) ^ . (DIST)^ . (RR)^(DIRl)^ 2RR-DIR1 -PM (1 ,K) 

a2 a2 ’ a^ 

2RR-DIR1.DIST ^ 2PM(1,K) -DIST (PM(2,K))^ (RR)^(DIR2)^ 


a^ a^ b^ b^ 

^ 2PM(2,K) •RR-DIR2 (PM(3,K) ) ^ - (RR)^(DIR3) ^ 


+ 2*RR-PM(3,K) ‘DIR3 _ ^ _ q 
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This equation is of the firm (ARR) (RR) + BRR*RR + CRR - 0 

(8.23) 


where ARR - iHSil! - 

sT ci 


(8.24) 


2 

(8.25) 


BRR = 2 ^(-PM(IrK)-DIST) *DIR l/a^+PM(3,K) *DIR2/b 
+ PM(3,K) •DIR3/c^J 
CRR = j^(PM(lRK))^+(DIST)^+2.0*PM(l,K) ‘DIS tJ /a^ 

“(PM(2,K) ) - (PM(3,K))^ /c^ -1 (8.26) 


Equations (8.24), (8.25), and (8.26) are evaluated by the 
program and (8.23) is solved to find RR. Substituting for 
the value of RR in Equations (8.14), 8.15), and (8.16), a 
point on the subreflector is obtained. 
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8.4. APPENDIX D 

DEVELOPMENT OF NORMALS ON A PLANE PANEL 
In the APRIN routine a certain nximber of perimeter 
points for each panel are read in. To determine a unit 
normal on each panel « the following procedure is applied: 

1) Any three perimeter points are used to form two vectors, 
as shown in Figure 8.2. 





Fig. 8.5. Formation of two vectors from three 

perimeter points 


where A = (Xj^-X2) i + 

A ^ 

B = (xj^-x^) i + (yi'ya^ i ^ 


2) The cross product operation is used to find a vector 

normal (N) to the plane defined by the vectors A and B: 



/N A A 

i j. k 


AAA 

i j k 

-► 

N = 

\ \ ^z 

= 

X1-X2 yi-y2 Z1-Z2 


®x ^ ®z 


X1-X3 y^-y3^ z^-23 


5= (yi~y2) * (Zi-Za^-^yi-ya) * ^ 

+ (Xj^-X2)'(Zj^-Z2)-(Xj^-X2) • j 

(Xj^-X2) * (yj^-y3) -(Xj^-Xj) (yj^-y2> k 
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3) The unit nomal N is computed by: N = — ^r— 

|S| 

4) If this normal on the surface of the panel has a 
negative x component, then the vector is inverted to 
yield a positive x component, since any normal vector 
on the surface of the reflector should be directed 
toward the origin of the reference system, i.e., along 
the positive x axis. (See Figure 2.5). 


8.5. APPENDIX E 


FILL ROUTINE FOR A VERTICALLY POLARIZED FEED 
The basis of this siibroutine can be found in (5). 
use it, the E- and H-plane patterns of the feed must be 
provided by the programmer in increments of 1°. 


f 



Fig. 8.6. Definition of angles 6 and e 
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In Figure 8.2« the angles used in FILL are shown. 
(From (5) .) 

-1 


(S s cos (sin B' cos 

where Q' and are angles in the feed system. 

and r, . tan'^ sine ^ °l!n»" 


(8.27) 


( 8 . 2 . 8 ) 


Figure 8.4 depicts the interpolation ellipse which is given 
by; 


U“ . V‘ _ 


where 


Hence : 


u * rcose, v = rsine# 

E - 90°' 180 


. . . ®e'-90 • %'-180 

r - Etot - — 5 ^ j 5 — jr 


(8.29) 

(8.30) 

(8.31) 

(8.32) 


The code of this subroutine is shown in Appendix F. In that 
code , PROJX = cos6 and PROJEX = sine . To insure vertical 
(i.e./0^) polarization, PROJEX is set equal to zero. That 
means u = rcose and v=0. Substitution for u euid v is 
Equation (8.32). 

Etot = V- 9° ^♦'.180 ^ . ^-»90 


^'= 180 • oos^c’’ 


COS e 


2 2 

where cos e = 1 - sin e . Since a 6 ' polarized feed is asso- 
ciated with the z component of a cartesian system P(3,I), 
and P(4,I) are given as: 

P(3,I) = Etot (along z) 

P(4,I) = 0.0 (along y) 


8.6. APPENDIX F 


LISTING OF THE CODE FOR REFLECTR 
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MAIN 

IMPLICIT MCALPM (A-M«0-2I 
neALPA MAJ0MCAI«MINC»IU8I*N0IIM 

COMPLCXPI6 KTOTia«4O0lftPI£LOVC4OQI«PI8U>2<A00l 

iNTficeR suMPCivsyiiPca 

CUMMON/PARAMS/AOPOMP«aeLLPtCCLLPtOl»f •PSl •PLNPta Ul •PLNUMNI Jl t 

• PE EOC 3 1 # ALPHA #08? A •GAMMA • JCLAM# AM# AOROR t •0CLLP8 #C8LLP2 • 

• PS12»r^S?2«POINri3l#HORM4JI#0UlirCI#NPNL#HPOINr#$aRPC8 
COMMON/APRPRM/NPTPPLtNPER iM 

COMMGN/COLOS/OCLT»XC« ANGINCtPMI J«4I#R0« XMIU ZMA»ZMN# VMX 
CONMON/COMTRL/NOPTO) •NLtST.lflPT# tCA80#l4.|ST <100 1 
C0MM0H/0IMLNS/r0<M«Z01M#VCTtZCT 
C OMMO M/EX TENT/Y M t N t VM AX • ZM I N • ZMAA 
C ONMO N/MA TH/P I • P 1 2 • PI 02 *0 TOR • R TOO 

COMMON/PATrRH/ErOT#AMIHORt3tOI«AMA/OR<0|#MINOR#MAJOR#NANGLE< 91 
OIMENStON P<0#27OOI .TPLOCTS I # 2PLO<70| •P0CRCT6 l•PR<2«0OO I 
DATA OONE/SHOOHE /#ML VL •HPARTS/O* 7/ 

DATA VL0«VH1 •2L0*ZH|/I •OO't'lO#-! *004l0 «1 •00410«*l *00410/ 

POOCX l•20•04UL0GI0<X| 

P0a<X»al0 •OPOLOGIOIXJ 

MAXPTS«2700 

CALL NPUTCP#NPAT) 

00 400 l«l •MPM. 

CALL APRTUR<P«1I 
PRIKT 777 

CALL OUAHTZ<P«MPERIMtI| 

PRINT 770 

lP<10iU«ll«EO*0) GO TO 00 
iSX«l 

IP (tOPT«EQ*ll lS0«-t 
CALL APRPLT<P*NPOINT#IS*J) 

PRINT 700 
00 CONTI NOE 

iP (YMLN.LT* VLOl VLO-YMIN 
IP <YMAX«GT«YH|| VHI«VMAX 
IP (ZMiN«LT*ZLO> 2LO-ZM1N 
IP (ZMAX«GT*ZHII 2H1«2MAX 
IP(1CASS.E0*|| NPER1MS4 
00 «» L«I»NPERIM 
PRU«LfMLVLI»P< l«L4NPOlNT| 

0$ PRC2*L4MLVL)»P<2 •L4NP0:NTI 
NLVL-MLVL4NPERIM4t 
PRI 1 • MLVL )«l *00440 
J JUM>0 

OO 200 K«l.HPAT 

CALL INT6R<P#MA3OR<K}«ANAjOR<Xl«ANINO0il«Xt*PIELOV#PIELOZI 
PRINT 779 
NANG«MANGLE(KI 
00 15 0 L«I«NANG 

£TOT< l•L♦l$UMl•ETOT< t *L4 14UMI 4P I ELOYI U 
ISO £ TtiT< 2 *L4 I GUM )«ETOT <2 #L4|SUM 1 4P IEL0Z4 Ll 

200 I 9UM»ISUM4NANG 
400 CONTINUE 
PRINT 7ol 

IP 4I0PT.U)*1I 60 TO 420 

IP ilOI(2 •U*E0*0l 60 TO 420 
YUIM«VHI<»YL0 


ZOIMbZHI-ZLO 

r cr»c vHi^YLO i/2 «o 

ZCT«( ZHI«>ZLOI/2«0 
CA4.L APRMAPiCRfNPNL •- II 
PRINT 762 
420 isgp-0 

00 770 l«ltNPAT 
NAN0 «NAn 6L£( 1 I 
F .^«AXVs-l • 0D«^40 
FMAXZ«-1 .00^40 

00 45 0 J« It NANG 

YFLOC Ji«COABS(EfOT( It J4ISUMII 
ZFLOC J) cCOABSCETOT C2t J4ISUM) I 
FMAXV«0NA XUFMAXYt VFLOC Jl I 
450 FMAXZ-ONAXI C FMAXZt ZFLOC J II 

1 SUM* 1 SUM4NANG 
0«AM1N04( 1 t 1 > 

FMTOfi* -60*000 
FMZ08« -60*00 0 
PMRMOB>-60«000 

P6R*4^ MAXZ4FMAXZ»FMAAV4FMAXY 
IF CFMAXY *GT •i•00-10l F M YOb«FOB ( F MAXV I 
IF (FMAXZ.OT«1*00-10) F'4Z0B> F OB ( FNAXZ I 
IF ( PM .GT*U 00-101 P6RM0B«P0B( PBR I 

PRINT 600 t MAJOR! II* AMAJUR! il* MINOR Cl I • C AM INOR C J • 1 I • J« I * 31 
600 FORMA TilHl *//y24X* 

• •TAMLE OF ELECTRIC FIELO STRENGTHS COBI *23X • 

• ///I «Xt •PRINCIPAL PLANE LF CUT IS • * A5 # • - ••FB.3** OEC* 

• //19X*«ANGLE •*A5.» FMnM».F8.3.« T0».F8.3** BY«*F6.3t* OEG» 
PRINT 666* MINOR!! I 

666 FORMAT!//13XtAS*4A. *00! Z/Z) • ,4X» *00! V/ZI * tAX. • OB ! Z/V I • *5X« 

• • OB! V/Y !• t5x • 'PgROa* •/ I 

00 70 0 it NANG 

PMERi K|sP6RMOB-lOO*ODO 

OBY «FMYOB -1 00*000 

OBZ «FMZO0 -100*000 

PmR«ZFLO( K)#ZFL0!K| 4YFLOCXI4YFLOIX I 

IF CYFLU! Xl.uT* 1 *00-151 OB V«FLd ! VFLO! K| I 

IF !ZFlO!K|.GT*1*0O-15I OBZ-FDB! ZFlO! K I I 

IF IP hR*G T* 1 • 00-20) P6EP!X|«P0j! P6RI 

IF (FMYOB.EQ *-60*OOOI ObYs-60«000 

IF !F MZi)B*EU*-oO *0001 0BZ«-60 *000 

OBZZ«OBZ-FMZOB 

OBYY> UBY-FMYOB 

OBZVaOBZ-FMYOB 

ORV2«08Y-FMZOO 

P6R0B«PmER!K I-PWNMOB 

PRiNf 690* OtOBZZ.OBYZtOBZYtCBYY «P6R0B 
690 FORMAT! 10X*F9*3*5F 1 1*51 
0«O^A M1NONC3 * 1 1 
VFLD! XI-03Y 
ZFlO! X)*06Z 
700 CONIINUE 

PRINT 750* FMAXZ *FMZOB*FMAXY*FMYOB 
FORMAT !// ISX *• MAX IMUM FILLO V AluES- •// 1 5X* 


750 


86 


1 • • 20LO6CMAACFietO-2l )a20C0C( • • 1FC16.7 • • • .OFF 12. 7// 1 tX • 

• • 20LMO«AX(FlftL.>-vn«?0L06C« • IPt. 1 9«7 • • .0PF12«7| 

pm NT 7S5* NPARTS 

I fS5 F0NMAT<//14X» 

• • lNTeAPO«.AflON NUM8CP US20 FOR iNTftORATION ••ISI 

PR! HI •MAJOR! II •AMAJQRCI I 

765 FORMAf!lHI«///20X»*PRlNCIPAL PLANE • ••A9»F7.3«« OEGREks* I 

CALL PLOTA(A4H NORMALIZED 2-COMPONENT OF SECONDARY PATTERN lOSI 

• •FMZ0R*2FL0*NANG«M1N0R|| I•AM1N0R!I• 111 
PRINT 76S«MAJOk! II .AMAJOR! I I 

CALL PL0TAC64^> NORMALIZED Y-COMPONENT OF SECONDARY PATTERN ! D6 1 

• •FNYOatYFLD«NANG»MlNOR! I l•AMlNORf 1 •! II 
PRINT 765»MA«ORf l|•AMAJORII) 

CALL PLOT4I64H NORMALIZED POMER PATTERN CORI 

• •PMRMOB«P«ER*NANG»MlNOR! II •AMINORf 1 « III 
770 CONTINUE 

IF f J0I!4.1I.CQ«II MR1TEI7.77SI 
IF !I0I!5«II«E0#0I STO** 

REMIND 7 
CALL PILIST 

775 FORMAT!* -1134*1 

776 FORMAT!/* FINISHED INPUT 

777 FORMAT!/* FINISHED APERTCR — » 

776 FORMAT!/* FINISHED dUANT IZ — — - • 

779 FORMAT! * finished INTGAT 

780 FORMAT! * 444 EXECUTED APRPLT 444 • 

781 FOf'IUI!/* PATTERN CCMPUTATIONS COMPlEI E 

782 FORMAT!/* 444 EXECUTED APHMAP 444 • 

STOP 

END 


09 


NPUT 

CUbAOuTiNt nPUT(P«NPA?I 
IMPLICIT RLALAd lA-M.O-Z) 

REA4.9S MAJUR(5)» M|r«OfMe»l»NORM 
C0MfA.£A»l6 £Tor(2t400) 

INTEGER &URFCUSURFC2 

COMROK/BlOCK 6/VC9L«2CBL •HFMABL* HFMI6L 
CUMM0N/^£c.0/EPI9ll«ET 19 I ) »NP «NT • XS«VS tZ ft 
CCMMON/COLU&/OELT«XC» ANGINCtPMC J. 4 I «R S« AMXtZMX.ZMN. VMX 
LUMM0N/C0NTRL/N0PTC3I tNLI&T tlOPT • iCASS* ILIST ( 1001 
COMMON/PARAMS/AORORF #aELLP»CELLP*DIST .PSI •PLNPNff Jl tPLNORMIJI • 

• FEEOI3) • ALPMA«atr A»GAMMA«XLAM» XX* A0R0R2«B£LLP2«CELLP2« 

• Pbl2*01ST2*POl NT C3I * NONMI 3 1 • SURFC 4 *NPNL*NP01 NT *SURFC2 
COMMON/PATTRN/ETOT* AM INORI 3* 0 I. AMAJOR I 5 I *M|NOR* MA JUR .NANCES I SI 
COMMUN/MA TH/PI *PI2 *PI 02 *OTiiM *MT CO 

DIMENSION P< S*27S0I* T ITLEI 40) 

DATA DONE/SHOONE / 

ICASS*0 
I OPT mu 

READ S* TITLE 
S FORMAT (lOAO) 

READ! I *1U I FEED* ALPHA *6ET A* 6AMMA* XL AM 
10 FORMA T( 7F 10. 4) 

IF (|CASS«N£*1) GO TO 3S 

READ! 1*20) SURFC2 .AONOR2 *0ELLP2 * CELLP2* 0 1 ST 2 • PS I 2* PO I NT • NORM 
20 FORMATi 24 *9X*SF 1 0.4/6F I U.4I 

JS HEA0CI.J7I SORFLl *NPNL. A0R0RF*8ELLP*CELI P*OIST*PSI*PLNPNT«PLNORM 

37 FORMA Til 4 */X*t2*9F i0.4/6F|0«4| 

IF I 1C ASS.Nt. 1 ) GO TO 40 

READ 11*391 i fPMCl * J) *l«l *3) * J«1 *41 

39 FORMAT (3F10*4| 

C EnO of MAIM REFLECTOR I nPUT LaTA 

CACL SUDPNTl Pi 
GO TO 43 

40 READ (1*41) I (PC I • J) *I>1 *3) • JBl *41 

41 FORMAT C3F10.4) 

C END OF SUB REFLECTOR InPuT OAT A 

43 4EAO(l*S0i AX •\rCOL*ZCBL.HFMABL*HMiBL 

SO FORMAT 4!»F 10.4) 

C FEED RA'JIATIGN PATTERN 

REAOC 1 *SSi EP 
READ! I *SS ) ET 
SS FORMA T( SF IS* SI 

READ CI*C»0l NOPT.NLIST 
1*0 Format C3ii«2x*isi 

IF INOPIC 11 .EQ.I .UR.NOPn 2) .EU. 1 ) REAO(1*70I i I L 1 ST C 1 ) * I • 1 * NL I ST | 
70 FORMAT (ICIS) 

1 SUM* U 
HP AJm 1 

rr M£AO( 1.501 MA JORCNPAT I • AMAJONCNPAT I .MINUMCNPAT I* ( AMiNOMC UNPAT I. 

• !• I • Jl 

5 0 FORMAT C AS .SA .FI 0 .4 • AS *SX • 3F 1 0 .4 I 

IF (MAJOR (NPA T) .EO.OONE) GO TO 55 

NANOLE(NI*AT M (AM|NOR( 2.NPAI )-AMlNUN( I •NPAT) I /AMI NOR ( J.M>AT|4 l.S 
li^ ( NANuCECNPAT » *GT • 7Si GU Tu 5S 
1 SUM* I SUM H*ANoLt ( NPA T | 

NPAT« NPAT »! 





90 



IF CNFAT.LT.6) 60 TO 77 

PRINT JJO 

STOP 

0:» PRINT 3J5 

STOP 

tttt IF CISUN«L£.400I 60 TO 95 

PRINT 340*1 SUN 
STOP 

VS NPAT«NPAr-l 

00 ^0 L« 1*1 SUM 

£T0TI 1 *L)«40«0O0*0«0O0l 
96 CTOTC 2*Ll«f 0*000*0«0001 

PRINT S76*TlTLe« ALA M*FE£0 .ALPHA *BeTA.6AMMA 
PRINT 577 *AC.yCBL.2CaL«HFMAaL*HFMiaL. NPNL 
IF CICASS.NE.il CO TO 160 
PRINT 576 

60 TO C 120*1 30.1 40*150* 160.161 I *SURFC2 
120 PRINT 579 *PO INT* NORM 

60 TO 179 

IJO PRINT 560 * A UR UR 2* SC LLP2 

60 TO 179 

140 PRINT 561 *AUR0M2 

60 TO 179 

ISO PRINT 562 *AOROR2 

60 TO 179 

160 PRINT 563* AOROR2.PS I 2 

60 TO 179 

161 PRINT 564*AOROR2*a£LLP2*C£LLP2*OI ST2 

179 PriNT 565 *CCPMCI .J)*l«l *31. J«1*4I 

PRINT 566 

160 60 TO 1220*230*240*250* 260*2701 .SURFCI 

220 PRINT 579 *PLNPNT *PLNOPM 

60 TO 300 

2 30 PRINT 560 • AORORF * BELL P 

60 To 300 

240 PRINT 561*AORORF 

60 TO 300 

250 PRINT 562*AORORF 

60 TO 300 

260 PRINT 563 *AURORF *PS1 

60 TO 300 

270 PRINT 564 *AUMORF *B4LLP.CeLLP*0ISr 

300 IF |nPnL.sR:. 1I 60 TO Jl 0 

1 OPT* 1 
NPNL* I 

310 PRINT 565 * i CPI I * Jl • 1* I • 31 • 3* 1 *41 

PRINT 567 
PRINT 569 
PR|7<r 600 *LP 
PRINT 567 
PRINT 566 
PRINT 6(^0 *er 
PRINT 400*NPAI 
OO 320 l.NPAT 

320 PRINT 500*MA JUNCM|*AMAJ0RCM|*6|N0P(M| *1 AM|N0N|AA*M|*AAal«3l 
330 format I •• ••••••••• ERROR-MOMt THAN 5 PATTERN •* 


OIU(UN^«iPAQE^ 
Of POOR QUAUTY. 
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33 & 

340 

400 

300 

376 


577 


576 

570 

560 

561 

562 

563 

564 


565 


56<» 

567 

586 

569 

600 


•CALwUCATIOHS 6£OUCSTEO 4666666466 
FOAMATi* 46»*464666 £RAOR-MOA£ THAN 75 AN6l.ES IN* • 

• one AATTEhN request 4646664666*) 

FORNATi* 4446666646 E‘iRUR ^ REQUESTED* • 1 5# * ANGLES TO BE * 4 
« calculated exceeds avail, storage 4666666666«) 

FORMAT!// 

« NUMBER OF PATTERN GROUPS REQUESTED. ••«••••••••«• S/I 

FORMA T< SA .AS* » «* .F |0.4*i0X«A5«* FROM* .Fi 0«4 • * TO*vFtO«4«* BT*« 
Fl0^4| 


FORMATllHl .////.tSXf* FAR FIELD RADIATION PATTERN CALCULATION */ 


//• ••10A6/* *tlOA6/* *«10A6/* *.10A6// 

• INPUT PARAMETERS* • •// 

• WAVELENGTH OF ELECTRIC FIELD* 

• LOCATION OF COORDINATE ORIGIN MRT FEED • .3F8.3 

/• FEED ROTATION ANGLESIALPhA.BETA.GANMAI JFB.3I 

FORMAT! 

• APERTURE PLANE LOCAT ION! XC) .F 7.2/ 

• SUB OiSH SHADOW CENTER COORDINATES IN APERT • PL.*.*«.2F7.2 

/• half major axis of sub OISH SHADOW tF7. 2/ 

• HALF MINOR AXIS OF SUB OISH SHADOW 

• NUMBER OF PANEL S IN REFLECTOR** •lA/ ) 

FORMAT (//* MAIN OISH DESCRIPTION AND ITS PARAMETERS* •/! 

FORMAT !• IT IS A PLANAR REFLECTOR • ./ 


A POINT ON THE REFLECTR SURFACE! X« Y«2 |« •••...*••«** .SFe. 3/ 


• COMPONENTS OF UNIT NORMAL TO SURFACE! Xt Y. X) .« • 
FORMAT (• IT IS AN ELLIPTICAL REFLECTOR 

/« MAJOR AXIS OF THE ELLIPTICAL REFLECTOR*** 

/• MINOR AXIS OF THE ELLIPTICAL REFLECTOR* ««**** * 

FORMAT!* IT IS A SHERICAL REFlECTR 

/• RADIUS OF REFLECTR 

FORMAT!* IT IS A PARABOLIC REFLECTOR 

/• FOCAL LENGTH OF THE REFLECTOR **•«*•**•.•...*• « 
FORMAT !*IT IS A PARABOLIC CYLINDRICAL REFLECTOR 

/• FOCAL LENGTH OF PARABOLIC CYLINDER*** 

* ANGLE OF ROTATION ABOUT X-AXIS ! PSI I «**•«*«<* • 
FORMAT!* IT IS A HYBER60LIC REFLECTOR 

/• MAJOR AXIS OF REFL* IN X DIRECTION* «*••*%..•• . 
/* AXIS OF REFLECTOR IN Y O IRECTlOH* .«.**.*...•• . 

/• AXIS OF reflector IN Z 0 IRECT ION**. ••..**««•. .. 

/• DISTANCE USED FOR TRANSLATION OF ORIG* OF AXES 

FORMAT! 


'*^6.31 

•• 

••F78.3 

>*4F6.3) 

• * 

FB«3) 

•• 

••FB.3I 

•# 

• tFa.a/ 
»«*FB«3I 
• • 

••F8.3 

>«*F8.3 

••F8.3 

>«*FB*31 


• 46646 PROGRAM IN SINGLE PANEL MODE 66666 * 

/« MiNlMUM*Y POINT ON THE REFLECTOR |X. Y »Z J •*••.«•....• .3F8 *3 

/• HAXIMUM*Y POINT ON THE REFLECTOR »3Fa. 3 

/« MAXIMUM-Z POINT ON THE REFLECTOR 


/• NIN1MUM*Z POINT ON THE REFLECTOR !X. Y *Z) «.«..«*.«*** *3F8.3 
> ) 

FORMAT !//* SUBOISH DESCRIPTION AND ITS PARAMETERS- */| 


FORMAT r//* PATTERN OF FEED IN ONE DEG INCREMENTS OFF-AXIS-«/| 
FORMAii* E-PLANE */l 
FORMAT!* H-PLANE •/) 

FORMAT C2X.5FI6.I0I 
Pl^OARCOSi-1 .0001 
PI2-P 16PI 
PlD2«0.56Pt 
OTOR«P1/180. 


RTOO« 180* /PI 

RETURN 

END 



1. 
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I 


APfIfUH 

St^BROUriNe APRTUR4P* 4C41.4. I 
IMPLlClf RLALPB (A-H»(>-2I 
REALAtt NNATC J|»NMAC»NORN 
INTEGeR SURPCl •SURPC2 
C OMMO N^AP RPR M/NP TPPL • NPLiU M 

COMMON/CAS5/SR«J>«£NC JI.AO* V0tZ0«V«i«RM«O*A02«V02*202»CR2C Jl 
C0MMCN/FLC0/EPI91 I •£? (911 ^NP • NT tXS •¥$ «2S 
COMMON/MATH/P1»P 12 tP I D2 •OTOH •RTU) 

CORROR/COLUS/OClT »XC« ANG1NC»PM( J«A|«RS« XMX« 2MX« 2MN. VM X 
COMMON/CONTRL/NOPrCJ) tMLl &T • I OP f • I CAS S» ILISK 100 I 
CPMMOK/PANAMS/AQRORP. BELLPtCELLPtOISr »PSI tPCNPHrC 11 •PLNOAM(J) • 
Pceuc Jit ALPHA. BE fAt GAMMA •XLAM»XXt ACM OR2*aXLLP2»C£LLP2« 

• P9l2t0IST2tP01Nl4 J) •NORM! Jl • SURFC 1 •NPNLtMPOiNT • SURPCt 
OIMCKSIUN AINVCltJI #B ( Jt 2 I • SB ( Jt 2 !• C( 11 • XI J|« A (!• J I •£ I ( J| • 

•PISt27S0l 

lPllCAS:»«eQ.II GO TO 10 
M«1 

OO 2 l«2tA 

|P(PCltMI-P(J»ll >3.2.2 
J M«1 

2 CONTI NUC 

XMX«P( 1 .Ml 
VMX«P ( 2.M I 
2MX-PI1.MI 

10 IP ( 1CALL*«.T« 1 1 GO TO 50 

ALPHAR«ALPHA AOTOR 
0£TAR*B£T AAOTOR 
GAMMAR«GAMMAPOTON 

A(t .1 l-OCO&C ALPHARIAOCOSCGAMMAR |-0S|N|ALPHAR|A0SIN(BETAR|P 

• OSINIGAMMANI 

A( 1 .2 l-OS INC ALPHARI AOCOSCGAMMANi^OCOSIALPHARI tOSlNIBt TAHI • 

• OllNCGAMMAMl 

A( I•J|■-OCOS4B£TAH)AOS1NCGAMMAR| 

AC2.1 l«*OSINC ALPHARIAOCUSIBETAR > 

AC2.2I* 0CO5I ALPHARI^OCOSIBETAHI 
AC2.il- OSINCBETARI 

AC J.l I-OCG5C ALPHARIAOSIMGAMMARI ♦OSINCALPHAR I AOSINCBLTARI# 

• ocoscgammari 

AC1.2 I -OS INC ALPHAR I AOS INC GAMM AH l-OCOSC ALPHAR | AOSINC BE TAR I • 

• OCOSCGAMMARI 

ACl.ll- OCOSCBETARIAOCOSCgAMMAHI 
00 40 1-1.1 

00 40 J« 1. J 
40 AINVC 1. JI-AC J* I I 

NPENI M-4 
NPTPPL-2000 

50 IE ClOPT«EO«0| CALL APR I N IP • I CALL I 
f MAX* 0*00 0 
r MiN-Pi 

55 PMIN-PI4P102 

PMAX-PI02 

dtt OO 05 1-1 .NPtRIM 

00 CO J- I . J 

bO X(J|-A1N94J.1|4PI1.1I#AIMVCJ.2|4PC2.I I4A1NVC J. II 4PC 1.1 I 

H-OSUMTCCXCI I4PEL0C1 I >442 ♦ IX 42 I 47 LLOC 21 I 4424C Xf J I 4PEE0C II 14 42 I 
PC 1. 1 l-OAHCUSCC XI 11 44 EEOC II I /HI 



I 

I 


I 
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»|N?Ml«nSlN(Pt !• I II 

IF (VINIhT«L?«I«O-I 0I »INrHT«l«O-10 
F4««l laFI-QAASlNU JUtl^FKtOCtl I/C HF»I Nf Hf II 
41 IF IFtUII.4T«TMAJI| rMAJC^Cl.ll 

IF CPtl*l l•Lf•TMINI rMIMBPCItll 

IF (PC/«II«4T.PICAI(| PMAA«PC/»|| 

IF CPC2«I l•Lr•PM|N| PHlNaP|4«ll 

44 CONTlNUt 

UCLP-PPAA-PMIN 
OALraf MAA-TMIN 

NPa040NfCOeLP«C>FLOAn NPrPPLI/OtLT|41» 0 
NP«CC Kf»-t I// 14241 
ANteINCaOeLP/COFLiiA • CHP) -2«4I 
IFCICAS4.CU.I I CAC.L F|NOXCCP«tl 
NfilJallCtf /<2«04AM4INCI*1 *0 
NT««4N?024| 

PM|NaPMlN-0.44AN6lNC 
PMAAaPMAX4a« 44AN4I NC 
rcr«l f MAX4fMlNI/I«0 
f Ml TC T- OF LUA f I Nf Ofi I 4AN6I NC 
iMAAarCf ♦OFLOAfCNfOll 4AN6INC 
00 9b ^1 vNl 
OJ 94 Ra I »NP 

PCI«NP4N|M4| J-l I4NPNC |a f M IN4 I J- II 4AN4 INC 
94 P42«NPCNI M4C J-l I 4NP4R laPMlNt 4R-I I4AN4INC 
NlNPaNl 4NP 
NPOINfaNPtNiMANr NP 
TMINalMlNAMrUU 
IMAAaf MAX 4PT00 
PMlNaPMlNAHrU) 

PMAAaPMAXFPTUO 

AN9lNCaANulNC4IITOO 

IFIlUlllf ILALLI«C0*II P*4lNf 1 97 t f MINt f MAA .PMI NtPNAJI* 

• ANOINCtNTNPtNPOINf 
107 FUNMAT4//* ILLUMINAriON 0A7A-V/ 

• • TMklA iLLOMINATiUN • »F9« J* • fU 

• •• F9.2/ 

• • PHI ILLUMINATION FHOM******** •••••••••••• ••••••F9«3t • TO 


79* J/ 

• INcMkMCNTAL ANCcf. 

• iNLHtFOMl TOTAL NUM0L9 OF 4CN£NATtO PATS / 

• total NUMbCN OF APCNTUNt PLANC POINTS** * ••IT) 


IF ( 40MFC I .Nl .SI uO lU 114 
C4P4laOLU44PSI4UTON) 

SNPxl a04l N4P4I40TUM) 

114 CALL F ILLPIP^NPOINTI 
00 1»0 0 l•l•NPOlNT 
4lNP«USINiP4/ •III 
LUSPaOCUSCPC /• I ) ) 
SIKI-OSIN4P4 I .1) I 
COSfaOCOSiPI I •! ) I 
004 I • I las INT 4CO^P 
00 4 *'f I I *41 Kt4 41 NP 
0114 Jt I |aCU4t 
01141 •/l■4FLi04l ) 

004 /• /l•»FLt04/l 


9k 


BBC Jt2l««^fEE0t Jl 
CALL MiJLT32(B«A»Bai 

60 TO (l20»l30»i40»IB0«l*0«l6il tSUBFCI 
120 ARaO.O 

BRaBU *1 lAPLHOflMCI l♦aC2•llORLNOIIlM2l4BI J« 1 1 OPLNOfliM 9 1 
CR»-CaCI»2I^PLNPNrC II lOPLHOfIMm 

• •CaC2»2I^PLNPNTC2l IOPLNORMC2I 

• -CaC3«2l#PLNPNTC3l l#PLM0RMC3l 
60 TO IBO 

I JO AR«B( 1 •! I 442/AORORP4424CBC29 I I4424BC J« l|B#2l/BeLLP4#2 

dR«-2«0*CBi I • lltBC 1»2I/AUR0RF4#24CBC2 •! I0BC2 *2 I4BC3 • 1IO0C3*2 11/ 

• 8ELLPt42l 

CR-BC I •2l442/AORURF4P24CB(2t2|4#24BC3«2l442 I/BELLP042-I «0 
60 TO lao 

140 AR-BCl 9114011 •! 141112 9 1I4BC29 1 14BC 3# i 1 4BC 3# II 

BR»-294C0C 191I4BC l92l4BC29ll4B(292l4BC39l |4B(392I) 

CM«aC 1921 4BC 192I4BC 29 21 4BC 29 2 14BC 3921 4BC 3921- AORORF4AOROIIF 
60 TO IdO 

160 AR^BC 2911401 29 114013911401 39 II 

BR«-2 *04 10129 114012 92 I4B1 39 1 1 4B1392|42904A0R0OF4aCt9 111 
CR >01 2921401 292l4ai392l40C392l44904AORORP4BCl •2l-4«04AOflOIIF440 
60 TO lOO 

160 AR«01394l40l39l l4CSPSI4CSPSl-2904B129ll40CJ9l|4CSPS14SNPSI 

« 401291 I4B1291I4SNPS14SNPS1 

0R«-2 904101391 I401392I4C6PS14C6PSI 

• -1 0139 1)401 29 2I4B1 29 114013 #21 I4CSPS14SMPSI 

9 4BC291 |481292|4SNPS146MPS|42904AORaRF4BC U 111 

CR-Ol 392l4B1392l4C6PS14C6PSI-2904B1292l4ai392l4CSPS14SNPSl 
9 4BC29 2I401 29 2l4SNPSI4SNP6|449 04AOflORF4CBCl92l-AORO#IFI 

60 TO lao 

161 AR^lOiU 1I442/A0R0RF442I-101291I442/BCLLP442I-CB1391I442/CCLLP442I 
OR— 2 9041 1011 91I4BC 1 9 21/ AORORF442 1-4 BCl9 1 1401 ST/AORORF442I-1BC 29 II 

9 401 2 92) /BLLLP442I-1013 91I4013 92I/CCLLP442II 

CR« 11 01 1921401 192I401ST40IST-29 04B1 1 9 21 40 fSTI /A0R0RF442I - 4 B1 2 02I 4 
9 01292)/OELLP442I-4043#2I4B1392)/C£LLP442)-|90 

60 TO 10 1 

100 IF llCASS9Ne9ll 60 TO 161 
IFIOAOSIAMI9LT9I 9OO-IOI R«CR/BR 

IF COAOSIARI.LT 91900-101 60 TO 165 
R«1-0R*06 6RT 10R4BR-4 9 04AR4CR 1 1/1 AR4AR I 
60 TO 106 

101 IF IOABSIANI9LT9I 9OO-6 I R«-CR/BR 
lFlUAaSlARl9LT9l 900-51 6U TU 105 
V«BH4BR-4 904AR4CR 
R«1-0R40S0RT 1 V) I/1AP4ARI 

165 CONTINUE 

X0«B1 1 91I4R-B11 921 

V0sO129l)4R-O12t2l 

20»B(39l|4R-B1392l 

IFl I96T.1 ) 60 TO 219 

AFl ICASS9E69 I ) 60 TO 219 

IFl ICALL96T 9I ) 60 TO 169 

IFl IOPT9£09l I 60 TO 190 

RI«06URT1 C AS401 I 921 |4424nrS4012 9 2ll4#2»i2»»B4 3921 I442I-I 9O 
THTMAX«0ATAN1-42S4B13 92I I/CAS4B1192II I 
f Nf A06»rMTNAX*2 9 54ANOINC4DIOM 
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KC*-< M ItOCUlbl THTAU6I»«( I tJII 
CUMST ■OAtt»|AC-«k I 
MQ TU 190 
109 ACaAS^CUNSI 
AA«1MN 

IF( ICALt.Wf. I ) liO TO AI9 
190 CALL riNOACtP.ai 
A AaiMN 

lPUUPI«t9«ll AC«AA 

219 au TO 42t0»230»t40»250«200*lol I tSUOFCl 

220 NHAfI lla^LNOAMl II 
NHAM2)aPLNUflMI2 I 
NHAfC JlaPLNOMMIJI 
uO TO 200 

2J0 r^ATI 1 l•-A0•OtLLPOOl/OSQllTCA0••0•oa.L^9•A♦| V0OOt#209#llOAOAOO099AI 
««HAr4 2l*-T09AUMCifir»<t/OSOMr C A099 2 9OCLU»#OA#CT0*«2*20OOt|0 

• AOflO»4f>A»AI 

NHAIC Jl ■•20* AUHOAF99. /OSiJOTjkC 0 0 2MCLLP90A#C¥009 2A20092|9 

• AOOOOF99AI 
S0Q TU 200 

240 NHAT4 1 !•- AO/AONOMf 
NHAT ( 2 I 0/ AUMOMO 
NHATl Jl *-20/ A OAOfl#^ 

S0O TO 200 

250 NHATC 1 •04AORU«ir/OOOMT 14 • 0 • AUOOfiro# t4T0002 tlOO 40 1 

OHATC 2I> -T0/DSUHfl4»04A0flOll044t 4VOOOt420442 I 

NHATl Jl« -20/OS0NTl4«04AOROHr44t4Y04O2#20442l 

GO TO 20e 

200 NMAG*0!MH 114 •04AOA0IIF OAuflORf ♦ 1 204C SPS 1 4 ONPSI - TOO SNP SI 4SNPSI I 042 

• ♦! V04ONPSI 4CSPSI-204CSPS14CSPSI 14421 
HHATl ll•2•0•AONORr/NMAL 

NHAT I 2 I*SNPS 141 204CSPS I- VOOSMPS I I/NMA6 
NHATl J|«C SPSl 41 Y04SNPSI-204CSPSI l/NMAG 
GO TO 200 

20 1 OeNaOSUAT I I I A0401 ST 14 42/ AIMOIIP444 I 4| ¥ 04Y0/BILLLP4O0I4 | 2 0 02 0/ 

• LeLLP440l| 

NHATl 1I«IA040IST i/I IAOIIOIIP442I40CNI 
NHATl 2)a-T0/l IOCLLPOOJ | OQtNl 
NHATl Jla-20/l ICLLLP442IOOLNI 
200 IP IICASo«NL.ll GU TU 209 
NHATl I |a- NHATl II 
NNATI2|a-NHATI2l 
NHATl J|a-NHAl IJI 

209 SCALAN«2«04|UI 1* 1 14NHATI 11401 2 • 1 1 4f«4A Tl 2 1 40 1 J 1 1 1 4 »»4A T I Jl I 
Oil 295 i «1 • J 

295 »N|Ll« lOlL • I I-SCALAM4NHATIL I I 
• T |BP| J* A I/M 
tPI 14 .1 l/M 

Cl I l•C05l 4i.U5P4L I 1- SINP4CPI 
CC2 laCUST 4>|NP4i T l4LUSP4tP| 

Cl J|a-SlNl4t II 
OO 40 0 N«I.J 
Ll INI aO.O 
OU 4410 N« U J 

400 t 1 INIaCI INI4AIN.MI4CIMI 

SCALAH-2. 0411 111 I4NNATI I 1421 I2I4NHAT I 21 4011 J I ONHATIJ I I 


y 6 


4>U SOO 

bOO LHCKIaSCALAA^HHATCRI-ftlClO 
IF C1CASS«N£.1I ru 5S0 
CA4.L CA&SAIPI 
PHASe«P12^(M^«IM^O)/IIL AM 
PI !•! I«V 
PI2«1 l«2 
PI J«1 l■CR2t2 ) 

P|4*l l■£H2IJI 
P|5*l )«PHASC 
6U TO 600 

660 V«V0»iXC-A0|A6MI2l/^l II 
2-206 lAC-AOIASMlJI/SHIl I 

JaOSORTII AC-A0IOIXC-A0I6IV-YOIO4 V-VOI ♦12- 20 I A 12-2011 
0 |P« 0 AaSlXC-AXl 

PHAS£*P12A|M#0^01FI/Ai.AMAPI6« II 

Pit# 1 l«Y 

PI2.I W 

PI J«1 I-2.RI2I 

PI4*I l-ftRI J) 

PI S«l l«PHA6£ 

600 CONTINUE 
RETURN 
END 


6UOPNT 

SUBROUTINE 6UBPNTIPI 
IMPLICIT MEALAOI A-N.0-21 
REAL# A NORM 
INTE6ER SURFCI.6URFC2 

COMMON/COLOS/OElT»XC* AN0|NC»PM| J.4I«M0* XMIU2MX»2MN» YMX 
COMIVN/PARAMS/AORCRE • BLLLP»CELLP *01 Of »PSI •PLNPN! I J I • PLNORMIJ I • 

• FEE0I2 !• alpha. BETA. OAMMA. XL AM. XX. AQROR 2.0ELLP2 .CELLPt t 

• PSI2.01&T2 .POI NTIJI .NORMIJI.OONFCI .NPNL. NPOINT . SORPC2 
OIMENSION P(6.2f50l 

00 JJ X«l .4 

RR-O&ORTIPMI I.AIAPMII .Ml APMI 2 .X ) A PMI2 .X I 4PMI J .X I APH| 3 .X I | 

OlRI»PM| I .Xl/RR 
0IR2«PM(2 .Xl/RR 
OIR3-PMI J.XI/RR 

ARN-OIRI AA2/I AUROAFA42I-IOIM2AA2/IBCLLPAA2I I-I0IR3A42/ICCLLPAA2I I 
BRR«2.0A| (-IPMI I .XI»OISTlAOIH|/IAORaRPA42| |♦|PM|2•R|A0IH2/ 

• 4 BE4uLPA42 II^IPMI J.XIA0IR3/ICELLPAA3II ) 

CNNaf CPM4 I .XI AA^ #01 ST A42 42.0 APM4 I .XlAQIST I/IAONORPAA2II- 

• 4PM42.X IAA2I/4BLIXPA42I-4 4PMI J.X| A42| /4CEIXPAA1I !•! .o 
NR«|-BHM4 0SORTIBRRAA2-4.0 AARMAcRRII/I AAR 4 AAR I 

•»4 I .X l«PM4 I.XI-RNAOIHI 
P42.X|«PM42.X|-RRA0|R2 
P I 3 .X I «PM4 3 * X l-HRAO I R3 
33 CONTINUE 
METURN 
END 


Vi Ji 
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suoAugriNE cassacpi 

INPLlCir REAL #6 CA-H«0-2I 
REALA a NHAT2C3I •MA6SR«NMA62«NLR«*RHATCJI 
INTEGER SURFCi»SURPC2 

COMMG N/PARAMS/AORORP • SELLP • CELLP *0 1ST «PS 1 tPLNPHT C J| • PLNORMI Jl • 
FEEOi 3)»ALPHAtBETA»6AMMA*XLAII«XX»AQRaR2»KLLP2»CELLP2« 
PSI2*01ST2*P0INT <3 l•NORMC 3 !• SURFCl»NPNL»NPOI NT* SURPC2 
COMMON/COLOS/OELT«XCtAHGlNCtPM(3*4l»RStXMX«2MX«ZKN«VNX 
COMMON/CASS/SR(3ltERC 31 1 X0« TO«ZO* Y»2 «RM .0 •XOE • V02 •Z0 2«ER2C3I 
COMMG N/MATH/Pl »PI2 «Pf 02 •OTOR«RTOO 
DIMENSION OCC3).EI2C3I»P(5»2750I*SR2(3I tCCZI 
MAuSR«OSORT(SR(l lASRU )ASR(2) ASR(2 |aSRI 31 ASH( 311 
00 5 N«lf3 

FIND DIRECTION COSINES 
OC(K|«SRC NI/MAOSR 
GO TO < I0«20»3U»40»S0«60) sSURFCE 
10 AA«0*0 

SO«NORM( I )40C( II♦NORM(2|AOC(2I4NORM|J I40CI3I 
CC-UO-POINT 4 1 1 J 4N0RM C 1 1 ♦€ YO-PO INTI 2 1 lANORMC 2l4 

• 4 ZO-POINT43I )4NORM43l 
GO TO 100 

20 AA-40C4I )4A2/AaR0R244294(0CC2l442/8EL( *2442 I440C4 3) 442/CELLP2442I 

ttBs2« 041 C X04OC4 IJ/AUR0R2442I44Y04OC42 )/0CLJLP2 442 I ♦ 

• C Z040CC3I/BEELP2442II 

CC«(X0442/AONOR2442|44V0442/BELLP2442)44Z0442/BELLP2442I 
GO TO 100 

3 0 AA«OC4l I4 0C4 1 1 40C42 I40C42 |40C 43 1 40C4 3 I 
BB«2. 041 X04OC4 II4V04DC42I4Z040C43II 
CC«X04X0#Y04V0^Z04ZO-( AOROR2I442 
GO TO too 

40 AAa0C42)4»240C(3|4A2 

BB«2* 04( V040C42 ltZ0 40C4 3|-2«0«AOROR24 0C4 1 II 
CC«V04YO^Z04ZO-4 4. 04AOROR24A2 I -4 4«04 AOROR2 4X0 1 
GO TO 100 

50 SNPSl 2«OSANIPS12 40TOR) 

CSPSI 2»OC0S(PS124OTOR| 

AA«(0CC3|40C43|4CSPSI24CSPS12 14 4 0C4 21 4QC 4 2) 4SNPSI 24SNPSI 21- 

• 4 2« 04OC4 2 )4DC 4 3) 4CSPSI 24 SNPSl 21 

B8-2.0 4Z0 40C4 3|4CSPSI24CSPSI242«04T040C4 2I4SNPSI 24SNPSI2- 

• 2.04 4 Y04OC43I 4Z04DC42I )4CSPSI24SNPS12-44«0 4AOROR240C4 II I 

CC« 4Z 0 442 144 CSP 512442)4 ( T0442J44 SNPSl 24 421-2* 04 4 Y04 Z04CSPS124 

• SNPSl 2 I- 4 4 .0 4 AOROR2 4AOROR2 )- 4 4 • 0 4 AOROR24XO I 
GO TO 100 

60 AA«40C4 II 442/AOROR2442|44 0C4 2I442/BELU»2442|4 40C4 3I44 2/CELLP2442I 

B8«2. 044 4 X04 0C4 1 ) /A4JNOR2442 I 4 40 1 ST240C4 1 I/AOR4JR2442)- 

• 4 Y04OC4 2 1/BELLP2 442I-4Z 04OC4 31 /CELLP2442I I 
CC«(X04XO/AOROR2 4 42 14 4 01 ST2442/AOROR2442 1-4 Y 0442 /BELLP2442I 

• -4Z0442/CELLP2442I-1* 0 

100 IF40ABS4AA).LT.1 .00-101 RM—CC/BB 

IF40AaS4AAI .LT.l * 00-1 01 GO TO 110 
V2aUB 4BB-4.04AA4CC 
IF4 V2 .LT.0.0 ) V2*0.0 
RM-s4-ad40SOHT4 V2I I/4AA4AAI 
1 10 CONT INUE 

X02«X04HM40c 4 1 I 
Y02«V 04RM 40C4 2) 


90 


Z02«Z04lUi«OCC Jl 

60 TO C120»l JOtl4O«l»O«l*O»l?0l»SUll^C2 
120 NHAT2C ll•HOftM( II 
NHAT2(2l«NORM(2l 
NHAT2C J)«NORMC3l 
60 TO 200 

130 NHAT2 Cl l■-A024BeLLP2•A2/OSOIITCA02•42•BCl.L^2A•44( T02442420J442I • 

• AOIU2I2444I 

NHAT2C2)a-Y02AAOH0R2442/OS0RTCX024424MLLI»Pt44A4C V024424 20244 21 • 

•AURQR244A I 

NHAT2 C 3»a-2024A0RUR24 42/0 SQR TC X024424aC4.LP24A I V 02 4 424 202442 1 4 
• A0R0R2444 I 
60 TO 200 

140 NHAT2C ll■-A02/AOROR2 
NHAT2 421 ■-V02/AOROR2 
NHAT2C 3)a-202/AUROR2 
60 TO 200 

150 *4HAT2C 1 1 >2 • 04AUROR2/0 SORT 44 • 04 ACROIU4 424V02 4424202442 I 
NHAT2 42la - V02/0S0RT4 4. 04AOROM24 424 T0244 24 2 024 42 1 

HHAXZ 43) ■ *202/OSORT 44 •04AuROR24424V024424202442 ) 

60 T(J 200 

1 40 NMAW»O50RrC4»U4AOROR24AOMOR24(2024CSPS124SHPS12- 

•T0245NP5I24SNP512I4424C V024&KPSI24CSPSI2-2024CSPS124CSPSI2I442) 
NNAT2 C 1 1-2.0 4AOROR2/NMA62 

NHA T2 t 21 aSRPSUO 4 2024C5P512-Y02 45NP5 12 )/NMA 62 
NHAT2 4 3|a C5PS 124 4 Y024 SMPSl 2-2 024C SPSI 21 /NMA62 
60 TO 200 

I 70 02R2aO&OR .'4 4 4 A02 401 572)44 2/A0R0R2444 ) 44 Y024 V02/atLLP2 444 ) 4 

• 4 2024 202/CU.4.P2444) I 

NH AT2 41 )a4X02401ST2 )/ 4 4 AOAOR2 442 I 40CN2I 
NHAT24 2|a-Y02/4 4 DCLLP24 42 1 40LN2 ) 

NHAT2 43 ) — 202/4 4 CELLP2442 )406N2 I 
200 5CALA2«2«044 0C4 1 )4NHA T24 I l♦OC42l4liHAT242l 40r 4 JI4NHAT2 43 ) I 

00 250 Lal»3 

250 5R2U.I aOC CL) -SC ALA24RHAT2 41.) 

£124N)aO.O 
00 300 Ra| «J 
300 L I24N laCR 4 R) /RM 
DO 350 A«1*J 

SCALA Ja2.04CL124 1 I4NNAT24 I ) 4L 12 4 2 ) 44HHAT 2 4 2 ) 4£ 124 3I4NHAT243) ) 

350 LM24R |aSCALA34RHAT24A|-C124A) 

1F40A0545R2C 1 ) ) .LT • 1 .OO-B ) SR24 lia|«0O-5 
Yavo2 44 AC-A02 )4 5P2 42I/5A2 4 1 ) 

2«2024CaC-A02|45N2C3I/SR24 1) 

0«0 SORT! 4 AL-A 02) 442 4 4 Y-Y02) 442*4 2-202 1442) 

HLIURN 

ENO 


99 


SUttMiVifINt APAINCP.ICALLI 
IMPLICIT RCALAS |A-H»U-2I 
HCALAa NOMM 
INTCCen &UMFCUSUMPC2 
COMMUN/APAPR M/NP f PPL • NPIRI M 

CailMON/PARAMS/AOAOPP*OeLLP*CLLLP»OIST tPSI •PLNPNTC Jl •PLNORMf J| • 

• rCEOI 31 .ALPHA. OCT A«6AMMA«XLAM*XC*AaRORa»6ELLP2*CeULP2» 

• PLl2»aiOT2.POINT( 31 .NORMC3 1 . SUMPCI .NPNL.NPOl NT .SUMPC2 
COMMON/FCCO/CPI9II.CT49I|.NP.NT.AS.TS*2S 
CONMON/CUNTRL/NOPTI 31 .NLI ST . ICPT . ICASS. ILIST C 100 I 
UIMLNSIUN PIS. 27901 

RCAOII.IOI NPCRIM.SURFCI .NPTPPL 
to FORMATI3I9I 

IF |NPCRIM»LC«2I CO TO 290 

IF INPCRIM.CT.40I CO TO 260 

IF ISURFCI.CT.6I C4i TO 270 
IF lNPfPPL.CT.2900» CO TO 270 

IF Cl NPLRIMASUMFCII.LC.OI CO TO 290 

RCAOC I .20 ) I CPI 1.3 I. 1*1.3). J-I.NPCRIMI 
20 FORMA TC3F I0.6) 

M*l 

00 2 1*2.NPER1M 
IF CPC 3.MI-PC 3.1 ) 13.2.2 
3 M*| 

2 CONTINUE 

AS*PC I.M) 

TS-PC2 .Ml 
ZS*PC 3.M) 

26 CO TO C30.40.90.90. 60. 611. SURFCI 

30 PLNUNMCI) *CPC2.I )-PC2 .2 ) I #CP 43. I I -PC 3. 3 ) I- 

• CPC 2. II-PC2.3II4CPC3.II-PC3.2II 
PLNLRMC2)*CPC3.I )-PC 3 .2 ) I 4CP C I . I )-P| 1.3)1- 

• CPC3. D-PC3. JII4CPCI .ll-PCI .2) ) 

PLNONMC3)«CPC I. I )-P| I . 2) ) #CP C 2. I )-P| 2.3)1- 

. CPCi.D-PCI •JI)4|P42.I)-PC2.2)I 

VMA^OSORTCPLNORMC I ) 4 424PLNOAMC 2)4424PLNOAMC 3)442) 

00 39 K*l .3 

PLNORMC 4- A) aPLNURMC 4-K)/«fllAC 

IFCPLNORMC I ) .LT .0.0) PLNURMI 4-K )*-PLNORMC 6»A) 

39 CONTINUE 
PLNPNTC ll*PC I .1) 

PLNPNTC2)*PI2.I ) 

PLNPnTC3)*PC 3.1) 

CO TO 100 

40 REA0CI.49) AUROMF.OCLLP 

49 FORMAT C2F 10. 3) 

CO TO 100 

90 REA0CI.99) AOMURF 

99 FORMATCFIO.J) 

CO TO too 

60 RCAOC 1.69) AORURF.PSI 

69 FORMAT C2FI0.3) 

CO TO 100 

61 RE^iOCI.70 1 AOMOMF.OELLP.CELLP.OIST 

70 FORMATC4FI0.3) 

too CONTINUE 


100 


199 

200 

2&0 

2*^2 

2C»0 

2o2 


270 

272 

J20 

3JO 

J90 

J»0 

JOO 

370 

401 

402 

403 

4 04 

4 Ob 

40o 


IFC 10 1C l»4CACLI«eO*OI HCTUIN 
PAINT 200* ICALL 

FOPIIArCM* •35«t*llAPL£Cr0fl PAMLL NUMMfl*«lAl 
OJ ru CJ20»JJ0«J40* J50t JAO»J70l«&UnPtl 
PAINT 232«ICALL 

FOAMATC///* 4944009400 INPUT CNAOA ON CAPC ONC FOA PANKL NUMMA* • 

• 14 •• EXCCuriON TCAMINATINO 4444444A44*} 

STOP 

PAINT 2o2»ICAU. 

FOAM4TC///* 0444444444 STQAAOC DOES NUT EXIST FOA NUAOEA OF • • 

• • PERI AE TEA POINTS SPECIF lEO - PANEL* tlA** 44444444S4t | 
STOP 

PAINT 2 72*ICALL 

FOAMATC///* 4444004444 AAAIAUA ILLUAINATION AEOUCST IS 2900* • 

• • MATS - PANEL*»I4*« ••44449449«//// | 

NPTPPL«2900 
60 TO 29 

PAINT 401 •PENPNT.PLNOAM.NPEAIM 
AETUAN 

PAINT 402 •AOAOAF«9ELLP«NPEAIM 
AETUIN 

PAINT 40J»AOAOMF.NPEAIA 
RETURN 

PAINT 404*AUR0MF .NPERlA 
RETURN 

PMINT 40&»AUACAF«PSI»KPEA1M 
RETURN 

PR 1 NT 4 Oo • AUHOAF • AELLP • CiAXP* O I S T • NPEA IN 


AETWAN 

FORMAT I/// IOXt*PANEL IS A PLANAR SURFACE**/// 

• A POINT ON TNE REFLECTOR SURFACE CX*T«2I*« 

/• COMPONENTS OF UNIT NORMAL TO SURFACE • tJFT.E 

/• NUMAER OF USER-SUPPL lEO E06E POINTS*** 

FURMATC///IOA**PA»iCL IS AN CLLIPTICAL SECTION**/// 

• MAJOR AXIS OF elliptical REFLECTOR *F7«2/ 

• MINOR AXIS OF ELLIPTICAL REFLECT OR ••••••••• ••F7*2/ 

• NUMEER of USER-SUPPLIEO E06E POI NTS** ••* • **I7| 

FORMAT l///IOX**PANEL IS A SPHERICAL SECTION*/// 

• RAUIUS OF REFLECTCR SPHERE**** ••**••* **F7*2/ 

• NUMAER OF USEH-SUPPL lEO E06E POI NTS* ******** * 

furmatc///ioa**panel is a paraoolic section**/// 

• FOCAL LEN6TH OF THE PARABOLA. «•••••••••••••••••••••* »F7 *2/ 

• NUMBER OF USER-SUPPLIED EDGE POINTS* ***** 

FOMMI TC///10X**PANLL IS SECTICN OF A PARABOLIC CTLINOER*/// 

• FOCAL CtNSTH OF THE PAH AOULA. ••••••••••••••••••••••* «F S* 3/ 

• FucAL LINE HuTAfICN FMUM V-AX IS C PS I I ••*••••••••**•* .FB* 3/ 

• NUMOLR of OSER-SUPPLIEO £06E PUI NTS. ••••*••••••••••*• 1 7 | 

fommatc///iox.*panel is a hvpemoolic section*/// 

• major AXl OF MEFL.IN X Ul REC T ION* .FB .3/ 

• AXIS OF MEFLECTOM in V O IMEC f I ON J/ 

• AXl^ UF REFLECTOR iN i DIRECTION 

• NUMBLM OF USLR-SUPPLIEO EU6E POINTS* 1 7 I 


ENO 


101 


FINOKC 


SUMUUflNC FlNOJlC(^»tl 
IMPLICIT ACAL^tf A-H«0-Z| 

CUMMON/COLClS/ueL f • AC f AN6I KC • PMC J *4 I tAS«KMII« /MJI« 2MN» V«A 
C0mttN/MATK/Pt«P12»PI02»0TCMI*flTO0 

COMMON/CONTML/NOPTCJ) •NLISf • lOPT. ICAftS* ILIAU 100 I 
OIMKNSION P<5»2700l #MCl»tl 
1P||CASS«NC.II GO TO 14 
M«l 

OO 2 1»2»4 

1PCPMC3»MI-PM|J»1I| J»2»2 
3 M-l 

2 CONTINue 

M^fMSuAT CPMCl •m; 442^PMC2 •.«I442^PMC J« HlOOtl-l .0 

TnTMAA«OATAN|-PM(3»M* <PM| UMM 

f HT AUG«THTMAA4J «04ANG INC 

XC«-M ^MI^OCUSC THf AUGI 

RETURN 

15 RSM«U50NTCi»''MA4liCl .21 |442 ♦ CTMX*MC2 «2 I I442^C 2MX^IM J. 21 |042l-l«5 

f H?MAA>OATANC-<2MA#BC J* 21 I/C AMX^BI 1«2III 
THT AUG«TNTMAX#k' .OPANG I NC40TUR 
2MN»-CR5M40CU5C TN r AUG I »BC1 •211 
RETURN 
END 


lOl 

/UNCTION lUl CINTENT*! fERI 
IMPClCir RtAL^B CA-H.O-/.I 
INTEGER 5UAFC1 *504/02 

COMMUN/CONfRL/NUPTC Jl »NC I SI • I OPT • ICASS. I LI ST C 1001 
GO TO CEO •30*40«;»0»4U !• INTENT 
20 IP CNOPTC 1I«LO«OI GO TO 00 

IP CNOPTC ll•E0•2l GO TO «l 

22 00 25 !■! •NLl ST 

IP C1LISTCII.E0«1TERI GO TO 91 
25 CONTINUE 



GO 

To 90 




30 

IP 

CNOPTC 21 .GT.OI 

GO 

TO 

91 


GO 

TO 90 




40 

IP 

CNOPTC 2I.EQ.0I 

GO 

TO 

90 


IP 

CNOPTC2I.EO«2I 

GO 

TO 

91 


uo 

TO 22 




50 

IP 

CNOPT lil .GE« 1 1 

oO 

TO 

91 


GU To 90 

IP CNOPTC il.EU.2l 
VO 101«0 

RET^N 
91 IU1«1 

RETURN 

ENO 


GO TO 91 


102 


SvIttMOJf INC tNf U<C^«MAJON*AMAJON« AII|NOII«F|CLDY»F|CU>CI 
IMPLICIT MCALMS CA-N«U-CI 
MCALt# MAJOMtMlNOA 
NCAL^A NORM 

CUMPLCX910 LTCMP.CCI •CC2«CVI •TS2«fST •Oil •OT| • 2IULO* VIOLO* 

• Z I vVl •PL4>2 •FLO Y*F ICLU2I 2001 •FICLOVC tool 
INTCCLN SUMFCa tSURFCt 
LJMMN/MATH/RI»PI2«PI02*0?0R»RIU0 

CCimN/CONTHL/NOPf|JI«N^ISr« lOPTt ICASS»ILIST( 1001 
CUMlV)N/P4KAMS/AuRORP • OCLLP tCCLLP *010? tPS I •RLNPNT C3 1 • PLNURM C J I • 

• PCLOC 3 l« ALPHA»OtTA»^«AMMA • ALAM* lU • ACMORttOCLLPt .CCLI-Rt • 

• PS 12 •0ISI2*PUINT (3l*Ni^RM4J)»SURFCl«NPNL»NRUlNf tSVlRPCt 
OIMCNSION AM1NORC3) •PCS*t?SOI 

OAT A HPHI •HTHTA/OMPMl •OhTMCTA/ 

SCN«VV9.0 

NPAMTSa? 

MPART«1 «0/NPAMTS 
2LAM-P12/ALAM 

CALL SETMCSCN.PC l•NPOINT♦l !•;») 

MAJOR 

OC6R-OLCAOTOR 

OLURvAMlNUMI 1 IPOTOft 

01cN«AM|MUN( J iPUTUR 

OSrOPR*AM|NURC2|POTOR ♦O ICR AO .9 

NTN«0 

0«IA.UR 

IF CHAJUN«Ne*MPM| I OO TO JAOO 
400 CO^«OCUSlOtoRI 
SINPbOSINCOCCR) 

L jsr»ocosioi 

SiHTaOSiNlol 
au TO J42S» 

J400 LQSPbOCOSIOI 
S INP^OSlNlOi 
CQST-OCUSIOCOJO 
;»lNl>0SlN40t0RI 
J429 NTH«NfN4| 

CTSP«CL;»T4S1 NP 
2A«ZLAM#COSr 
V K«a. AM4S INP4SINT 
I OLD* I 
INCm«2 

FL0T« ( 0«U .O.O I 
FLOZ^I 0*0 •0« 01 
VULO-SCN 
VI *4 0*0«0«0l 
2l«40.0.0«0l 
J490 COKflNUC 

IF 4PI I •lOCJ I «Nt •PI I • INCAI ) OL fO 4000 
2«P42 • lULUl 
I. R V ^ I J • I OLD I 
fc.R2«P| 4.IULU) 

PH«PI9 • lOLOl 
J2»4PC 2»lNt BI-ZI4MPARI 
OLMT« 4P4 J* INLM I-LRTI4RMANT 
OLRZ-4P44 • lNCal-LR/l4NPA^* 
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4> C P I S» 1 1 NC ■ ) - PH I • IIP AA f 

CriLI«*aCOCJU»COCMPLIUO« OOOtiKA^-PMI I 

Cil«tN2»COSPACTCMP 

CVI«(CArASlNf AfJI^PCfSPlACTCIIP 

Tsr«c o«o»o*oi 

rs^«io.o»o«o) 

Oil J700 N«l.NPAAft 

CAV^fcRV^OLMV 

€RZ»tRZ^tMl 

»»h«pm#opn 

CrCMPaCOeXPCOCMPLACO. 000*2AAi-PHl I 
Ci2«CMZRCOSPACriMP 
C SI Nf ♦CII^ACf SPI AC? CMP 

TSZ«TSZ#CZI^CZ2 
TSV^f Sr>CTl ♦CV2 
Cil«C22 
CVI«CV2 
J700 CUAflNUC 

♦TS2AIO«SAOZI 
r i«Ti ♦rs? •lo.SAuzi 
J900 l4ii.O« lui.04 1 
iNfcV* INIs A1 
60 TO JAftO 
4000 ^OHTlNUt 

TNAVaPC !• lOtOI 
IP C VOLil.CO.SCNI 60 TO 4400 
4200 02l«( 2I-2IOCOIAAPAIIT 

07I-C 7 1-y 10LO|4flPAMT 
«>V«4 YNLO- VUL0I4RPAIIT 
cr AMP«Ci>tAP(OCIiPLAC 0«000» TA4YOLOI I 
L2I «21 uL 4>4CTCMP 
CTI«T lULOACrtMP 
Tsy*io«o«o.oi 
r S2«t o*o«o«oi 

OO 43yU N«l«NPAAfS 
yULO«VC4.0#0V 
ZiOLD"ZiOiM^OIi 
y luLOay lOi^yoYl 

CrCMP«K.4>4.AP(OCMP4.ACO«OOJ»yii4yULOI I 
C22«2 lOLUACTEMP 
crz»f iOLO^crtmp 
f S2«IS2 4C2I»CZ2 
f sy«T sY4Cy I »CT2 
C21-C22 
cvt^ V2 
4J00 CUNfiNUi 

fLOZ^t Lo/^ rs24C «s40yi 
PLOr«iPLUV ♦TST 4| • 640 Vi 
4400 CUN7I NUC 

y«iLi>« y^4 

2 lcM.Ua/1 
VI iiu>«yi 
V 1 «I0 •«!« tf • 01 
21 «cy .OtO .01 

IF CPC A •! NL «l .Ni .SANI 6U fO JvOO 


5000 


5 

10 


100 


FUL0V(NTM|«FLDY 

F1ELD2(NTN)«FL07. 

D>D^D1CR 

IF (D*6T«0STOPR) 60 TO 5000 

IF (MAJ0R«5Q.HPHI> 60 TO 400 

GO TO 3400 

CONTINUE 

RE TURN 

END 


fillp 

SUBROUTINE F ILI.P (Pf NPTI 
IMPLICIT REAL^B (A-N.O-21 

C0MN0N/FEE0/EP(91i •ET49II tNPtNT •XStYStZS 
CONMON/MATH/P|.P12«P1D2»OTOR«RTOO 

CORMON/CONTRL/NOPTI J I tNL 1ST • tOPT • IC4SS* IL 1ST 1100 1 
DIMENSION P<5*NPT) 

00 too Isi.NPT 

PROJAsOSlNiPC i .1 n«0C0SCP<2» 1) I 
PR0JEX«O«0D0 

IF 4 0ABS<P(2 f ll-Pl I.GT*! .OD-S) 

• PROJEX»OSINC OATANlOCOSCPi 1 • 1 > »/DSlN(P(l » 1 n/0SlN4P42« XI I I I 
ANGLX«OARCOS(OABS( PROUXI l•RTOO 
LO«4NGLX#i*000 
IH1«L0^1 

PPFLO«4AN6LX-orLOATCLO->l J ) # 4 f P 4 1 Hf)-EP4 L O )) »EP CLO I 
TPFLO*4ANGLX-DFLOATCLO-1 MP4ET4 IHII-ETC LO) ) «>ET 4L0I 
SI 1C 2« PRO JEX • PRO JEX 
COS£2«l«0OO-SlNE2 

PI 3* Z I•PPFL0•TPFL0/4 0SQRT(TPFLD«TPFL0•C0S£2♦ 

. PPFLO#PPf L0ASINE2n 

P4 4. I l•0•000 
PCS* X )«0 .000 
CONTINUE 
RE TURN 
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OUANfl 

;»UI»RUUTINE OUANrZ(P»NPCiUM«lCALLl 
IMPLICIT HeALPM CA-M»0-ZI 
MCALPd NORM 
iNTtfuLR &UHFC1 •^UAPC2 

C UMMU N/dL OCK G/ TC BL • ZC BL • HP MAUL • HP Ml IIL 
CUMMUi^OI MLM^/Vd|M«ZOIM.rCT» ZCf 
CUMMUN/eXTLMr/yMINtVMAA.ZMlNvZMAA 

CCMMON/CUHTRL/NUP743l«NLI$f • lOPT • ICASS* ILISTI 1001 
CJMMUN/F£C4>/EPC VI »NP •KT .XStVS *ZS 

COMMUM/PARAMS/AuRURF«BeLLP.CELLP»UlST »PSI •P».MPNT( 31 •PLMORMI3I • 

• FEeO(3l«ALPHA«BCrAtGAMMA»XLAM«XX»AOROR2.aCLLP2«CELLPt. 

• P:»l2tOISf 2»POINT|3l .NORMIZI •SURFCI •NPNL.NPOINT .SURPCZ 

0 1 MENS I ON PC5»NP4ilNTI«PINT(SltPOLOf Bl tPBLKC 5l*PRB<BtAl» •ZC2»I0II 
IFC ICASS.EQ«II NPbRIM«4 

NBARSbNP-2 
YMIK«I .OO^IO 
VBAX«*I.0U^10 
ZM|»<^I •OO^IO 
ZMAX«- I • 00^1 0 
NOS«2 #NBARS 

CALL SETMC I •0O^20«ZtNOSI 

00 20 I>I•NPERIM 

IF (PCI»I l•GT•yMAXl YMAX«P4I.1| 

IF 4PC I tl I.LT^YMINI VM|N«P4l«n 
IF (P (2* I )«GT •ZMAXI ZMAXsP(2«ll 
IF |PC2»I I •LT.ZMINI ZMIN-P42»ll 
20 CALL MOVLMIPC !• I I »PRS4 I tl I •SI 
YOIM- YMAA-YMIN 
VCr«4 YMAA#VMlNl/2« 

ZOII^ZMAX-ZMIN 
ZCraC 2MAX#ZMlKl/2 • 

GRID* 4 YMAX-YMINl/4OfLOAnNBARSI-0«6| 

GRIOLObYM 1N#GRIO/S *000 
CM IOHl«VMAX-UMi J/S.OJO 
31 IBGN«NPCH IM#I 

NOEX«NPLRI M 
JO 100 l«lliCN«NPOlNT 
IF (P( 1 • I I.GT .VMAXl GO TO OB 
IF IP ( I .11 .lt.vminI gu to VB 
NGRIO«(P( !• 1 l-GRIULOI/GRIOfO«S 
P(l •! I«GR1JLOH>FLUAT4NGR10 I*gR 10 
CALL MOVEMIPI Ul I .PI I •I-NOLXI »5I 
GO TO 100 
Va NOLXaNOLX4| 

100 CONflNUL 

NP01NT«NP4ilNf-NO£X 
CALL PISUNICP.S.NPUINTI 
IF ClUPT.LO.il GO TO A22 
CALL MOVEMlPMSCl .11 .PRS Cl .NPERlMtll.BI 
X OL X« 2 
Y2«PNStl • 1 I 

2 2«PN SC2 • II 
Y I^PNUI 1 .M>LXI 

1 t«PNSC2 .xucXl 

IF COABSi Y I- Y2I .L T. I • 0O>SI GO TO 400 
SlOPL«I2I -22 I/4Y 1-T2I 


200 


106 


IF IVI-Vtfl MM094tJO$dJO 
idO VNl«V^ 

V4.0«V I 
UU fU i40 
dSO VHl«VI 

d^O lNa&A«l VLU-^lU4.0|/(iMll>»l«0 

VO«<iN |IK.U♦il^ tUA T 1 1 NOt Jl» II 06H l4> 

IF (Vg«Gf.VH|| GU TU 400 
iLilAil*! 

ZC£*:»4.ilFlL •VU#« 

IF 1211 tlNUk4|*Lf«l«0D4l0l l4.UAO«l 
21 ILOAU«lNU«.«l«2Aft 
GO TU WtMl 
400 V2«V1 

22«2I 

A0AA«IU>t:A4l 

IF IA4ICA«LI.*NF1N|M4|I GO fO 200 
OJ 420 l•l•N•AHS 

IF 1(21 I •Il4il2tl I I.Cf •! .OO^IO) GOTO 1004 
IF I2I2*I l*2( (•III 4IOt420»420 
410 22«2(2«ll 

242«l 1«2( 1 •! ) 

211 tl 1«21 
420 CilAflNUC 
GO ru 444 

422 HF MAL A« VO I *%/2 •OOO 

HFMIl )I^OIM/2«000 
00 4J0 1* I • NOAMS 

V«gMI 0LO4OFLGAI ( 1-1 I4GII10 
VA« UUOO-C C V-VCT|/HFMAtM|442 
IF ( VJ .LI «0«0 I VJ«0«0 
II MI4A4OS0NI ( VJI 
21 1 • I l•-22♦2CT 
2(2*1 l« 2242CT 
4 JO CONtiNUi. 

444 L«0 

N«1 

CALL SLIM! 0«0»POLAtSI 
VO«P( I •! I 

lOfcA«0INT 4 ( VU-GNIULUl /GNI 0» 1 .001 I 
00 OOO 1 K UNPUINT 
If (PI 1 *1 I •LO«VO| GU 10 4o0 
IF (L «gT *2 I «#0 10 4 fO 

N«N-L 
4/0 L«0 

VO«P( 1 • II 

I0FA«0LNI ( ( TO-GIIIULUI/GMI04* *00 1 I 
POLAI 1 (•PI I • I I 
P0LA(2I*P|2*I I 
U PT«- 1*0 

IF IP(2«ll.eo.2( 1. IOlAl«OM.P(2tl l•tO•2(2•IOtX|| VIST* 0.0 
IF (P(2*l I .v^T •21 1 •lOCAl •AN0.PI2 • 1 I •LI •2I2»ICH;A| I Il4f»l«0 
rLSfUL«»0 MAOL 4HF MAUL 4Hf M I BL 4MrM10L 


400 


t>ol 

0 

l:> 

ro i 

7 04 


71 I 


7J0 

7^t> 

OO J 

VOO 


Vt>0 


O J I 

10 0'^ 
\ 0 JO 


~HFMAtiL#HFMAttL#CP|2 •l)-7CBC»*(Pi^ .ll-ZCaLl 
-HFM lUL4HFMl8L»iP< 1*1 )-VC6LI*(P( 1*1 |-VC»L> 

IK lTt:»T) 70t*&0l«501 

iF (T tSTaL«Lt«0.0> UJ TO 510 

CAUL HUVt MiPBLK .PC t «NI t5) 

oJ Tu :>15 

call mUVLM(P< 1 • i ) dP(1 *N) *5» 

N*N* I 
L »L4 I 

IF iTLST.LO«0*0l OO TU 000 
IF iL «LO*OI OO TO 600 

IF iTF^i^TEsro) 7o4«aoo«aoo 

lALX iNTPLiPCU) «P( I *i ) tPlKT .Zil • 10£X> I 

NL«0*C 

iF C T LST. LT. 0*01 OO TU 7t I 
CALL MUVLM* I 1 I .Pi I *N| 


NCMG* I 

CALL HOVb Mi PI NT .PbLA.^I 
T £ST 6 L «HF MA8L 4HF MA6L AMFM 16L4HFH lOL 

-HFMAULAHFflASL# iP I NT 1 2 >-ZC8L » • i P INT ( 2 ) -ZCOL ) 
- HF Ml BL4HFM I BLA (P| NT< I )^YCBL) • (PINTU )-VCBL> 
iF ( Tl STUL«LL *0 *0 ) GO TO 720 
CALL MOVt MiPBLK #P4 I .N-NCHGI #51 
CO TO 725 

CALL N.LVtMiPiNT*F4l •N-NCMG) .5) 

N»Ntl 


C «L ♦! 

CALL i40VL MiP ( 1 • i I •PULL) *51 
T tSTO - Tt:>T 


C uNIl NUt 
NPO IN I *N- I 


NLUCc l>#NPLMIM 


CALL MCVbMiPR5*PU «N| .NLUC) 

IFl 10 U4# i ) .tO.l ) liKI ILi 7 NPOINT.ILAUL 

1 Fi iC»i 14 • I ) «LO. 1 ) «M1 TL(7.9li2i i ( P i I • J I • I > 1 • 5 I • J> I • NPU 1 N T I 
IF t lu l 4 i • IcALlI .LU* 0) MlToNN 

PRINT ^50 • t MIN. YMAA. CMIN.ZMAX, uN lOLU.GNlOHl .GAIO^NBARS.NPOI nT 
KJ«MAl(//» 0UANTI2INO OATA-*// 

• POiNT PATTLRN tXTLNTS ON APtRTURL PL A N£ . • • . . • • T M I N* • «F 7* 2/ 

• • • •• • • .YMAX« • .F 7. 2/ 

.ZMIN-* ,F7.2/ 
2MAX«« tF7, 2/ 

• oHiJ NANOtb from*. *.F0#3«* To 

• , F a. J/ 


:>PAclNv. OLTWLLN y^RlO aANl. IS 


tmcrlforl NUMULR of s*K (0 bars lo/ 

NOMOLH OF POINTS SOPPL i L O TO RADPAT.^ *.17 ) 


) Oi4MA 1 42 1 bi 


) OMMA1 4 SO lo « t J ) 


K U T OR N 
PRINT kOOti 
F ORMA 1 4// / 


♦♦ I AMKAY FOR POINT RtSiOENCt NOT FIlLLU CORRECTLY 
//lOA.*- STOP LACcOTION -• » 


T uP 
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8.7. APPENDIX G 


OUTPUT FOR TEST CASES A AND B 


ORIGINAL PAG£ IS 
OF POOR QUALITY 
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